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AERODYNAMIC LIFTING SURFACE THEORY 
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ABSTRACT 


This document Is a description of and users manual for 
a USA FORTRAN IV computer: program which evaluates 
spanw se and chordwise loading distributions^ lift 
coefficient^ pitching moment ^coeff Ictent/ and other 
stabil Ity derivatives for thin -wings In linearized^ 
steady^ subsonic flow* The program is based on a 
kernel function method lifting surface theory and Is 
applicable to a large class of planforms Including 
asymmetrical ones and ones with mixed straight and 
curved edges. 
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1 INTRQDUCtmN 


This document Is a description of and users manual for 
a USA FORTRAN IV computer program which evaluates 
spanwtse and chordwise loading distributions/ lift 
coefficient/ pitching moment coefficient/ and other 
stability derivatives for thin wings In linearized/ 
steady/ subsonic flow* The program Is based on a 
kernel function method lifting surface theory and is 
applicable to a large class of planforms Including 
asymmetrical ones and ones with mixed straight and 
curved edges. This program is used In conjunction with 
other/ separately documented programs which (1) set up 
the geometry (ref. 1)/ (2) set up the boundary 
conditions (ref. 2 ), (3) determine the aerodynamic 
Infiuence matrix (ref. 5)/ and (4) solve a system of 
linear equations Cref. 4). The aerodynamic theory Is 
described In ref, 5. 


(luestlons concerning either this document or the 
computer program or the associated computer programs 
should be directed to 

R. T. Medan 
Mall Stop 247-1 
Ames Research Center 
Moffett Field/ 

Ca. 94035 
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2 PROGRAM DESCRIPTION 

The normal loads program evaluates spanwise 
distributions of loading# lift coefficient# pitching 
moment# and center of pressure; the chordwise lifting 
pressure distribution; lift# Induced drag# rolling 
moment# root bending moments#, and pitching moment 
coefficients; and vortex drag factor# spanwise centers 
of pressure# and the lift on each half of the wing. 

The program uses what Is called command format 
programming. With this type of program the user 
himself controls the program flow calculating Just what 
he needs and In the order that he wants to calculate 
It. For the most part the card data Is entered In the 
order that the user chooses and .always In a standard 
format (8F10.0 or 16IS). This makes the program very 
easy to use. A description of the available commands 
is given In section 5. 

The program has commands allowing the user to select 
the spanwise and/or chordwise stations at which he 
wishes to evaluate any of the various distributions. 
The user can choose sets of default stations (no 
action required by the user) or he can enter various 
types of parameters causing the program to 
automatically compute sets of stations# or he can enter 
tables of stations manually. 

Although some data Is Input from cards# the majority of 
the data that the program requires Is stored on two 
disk or tape files: the geometry file and the solution 
f I le. 

The geometry file contains all of the geometrical data 
such as the local chord distribution# the aspect ratio# 
and the longitudinal reference length. This file Is 
created by the geometry program (ref. 1). 

The other file Is the solution file# which contains the 
coefficients in the expression for the lifting pressure 
coefficient. This file Is created by the equation 
solving program (ref. 4). Several sets of coefficients 
may be contained on this file. For example the file 
may contain sets of coefficients corresponding to (1) 
angles of attack# (2) pitching about some axis (for 
computing quasi-steady pitching derivatives)# (3) baste 
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camber and twist distribution^ and (4) rolling (for 
computing quasi -steady rolling derivatives). 

Each set of coefficients Is termed a case. The program 
works with linear combinations of cases called 
combinations (because Induced drag« for example/ Is 
nonlinear). The program uses a matrix of weights to 
determine the combinations from the cases. The columns 
of the matrix correspond to the combinations while the 
rows correspond to the cases. If the default option Is 
chosen by the user/ the weight matrix Is computed 
automatically and Is an Identity matrix. Otherwise the 
user enters the weight matrix. 
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3 


DEFINITION EQEFFJCIgNTS 


This section presents the definitions and expressions 
for the various coefficients calculated by the program. 
Fl|ure 1 shows how the various geometrical quantities 
referred to below are defined for a yawed wing. 

PRESSURF COPPCiriPMT 


ACp 


• •'u 

NMAX KK 


• 2b/c 


BNK(N^K)*HN(N,0)*SIN(K*e) 


N-1 K-1 


where 

b * the effective span (see fig. 1), 

c - the local chord (a function of eta), 

eta • the nondimenslonal spanwise variable such 

that -lietaji^l^ 

0 • -C0S"lx, 


X “ the nondimenslonal local chordwise 

variable such that -llxi.1^ 

6 • ♦COS“^eta# 


2./PI*COT(0/2.) 

or 

HN(N,^)- 2./P|*SINC(N-l)*$) for N>1, 
and 

PI - 5.141592... 
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SECTIONAL HFT CaEFF_LC.I£KT 


CL 


♦ 1 



ACp*dx 


KK 


■ 2b/c 



K-1 


SIN(K*8)*B(IC) 


where 

B(K) - BNK(l,K)+BNK(2,K)/2. 

(B Is denoted by STORl In the program. ) 
NONDIMEN&LQNAL 01 ACULATLQJj 
GAMMA - CL*c/(2*b) 

KK 

- SIN(K*9)*B(K> 

K-1 

NQRMAUIEP-LQAailLG 
CLCCLC - CL*c/(CCL*CAVG> 
where 

CCL • overall lift coefficient 


and 

CAVG ■ average chord 


5 



NORMAL LOADS PROGRAM 


JlCLNQtMEMS±QNAL PITCHIWG MOMENT 
QMOMNT - 2*cm(l/4)*c/b 


KK 

- *5* ^ SIN(K*9)*(-BNK(2,K)+BNK(3^K)) 
K*1 


where 

cm(l/4)» 




♦ 1 

ACp»(x^.5)»dx 


-1 


CENTER _QF PRPSSURg Pi S litlAUTtOM 

This center of pressure is normalized by the local 
chord and Is measured from the local 1/4 chord. 

COP • -cm(l/4)/CL 


- -,25*QMOMNT/GAMMA 
(LV.E&ALL LIFT CQEFPICIEMT 


6CL 


• 1 


.nff. 


ACp*d$ 


WING 


• AR/(4.*BRATIO**2) 


♦1 ♦! 



e* Acp*dx*d(eta) 


1 -1 


* PI*AR/(2.*BRATIO**2)*B(l> 

where 

S « wing reference area/ 

AR ■ aspect ratio 

- 4.*BREF**2/S/ 
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6REF ■ lateral (spanwlse) reference length 
(usual 1 y b/2)^ 

and 

BRATIO - BREF/(b/2) 


(MUCIQ .QBAfi 

The drag ts calculated using a Trefftz plane analysis. 


GDI 


KK 


PI*AR/4.* 



K*B(K)**2 


K-1 


VORTEX DRA&_F^T.QR 

The vortex drag factor ts the reciprocal of the 
ideal Oswald drag efficiency factor. 

VORD • PI*AR*CDI/CCL**2 

- 1 ♦ 2*(B(2)/B(1))**2 

♦ 3*(BC3)/Ba))**2 

♦ 4*(B(4)/B(1>>**2 ♦ ... 


ROLLING MOMENT 


CMR 


l./(S*2.*BREF) 



ACp*Y*dS 


WING 


» PI*AR/C8*BRATIO**3)*B(2) 

where 

Y ■ dimensional spanwise coordinate 

■ b*eta/2. 
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RIfiHT 

CMBP 


LEEl 

CMBM 


AQOT BewntMft 

• the nondtmenslonal moment required to 
hold the portion of the wing on eta>0 
I n equ 11 1 br I urn. 


■ 1 ./ 


(S*BREF) 



right half 


ACp*Y*dS 


- AR/(4*BRATI0**3)*(^P|*B(2)/2 + 


KK 


^ B(K)*(SIN((K-2)*P|/2)/(K-2) 

K-1 

Kt2 


-SIM((K+2)*PI/2)/(K+2))) 

OaT BENDIMft Mnur^tT 


the nondimenstonal moment required 
to hold the portion of the wing on 
eta^O In equilibrium* 


■ 1./(S*BREF) 



ACp*4*dS 


left half 


• AR/(4*BRATI0**3)*(-P|*B(2)/2 ♦ 


KK 

V" 

/ B<K)*<SiN((K-2)*P|/2)/CK-2) 


K-X 

K+2 


-SIN((K+2)*PI/2)/(K+2))> 
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L I PT J[1M_RLGHT S I GE _Q£^MtAG 


CLP 


■l./S 



^Cp*dS 


right half 


- AR/<2*BRATIO**2)*(Pl*B(l)/2. + 

KK 

B(K)*(SlN((K-l)*Pl/2.)/(K-l) • 

K“2 


SIN((K+l)*PI/2.)/(K+l))) 


L LFT ■Q^^ LEf_t 


CLM 


- 1/S 



ACp*dS 


left half 

■ AR/(2*BRATIO**2)*(PI*B(l)/2. 


KK 

^ B(K)*(SIN((K«l)*PI/2.)/(K-l) - 
K-2 


SlN((K+l)*PI/2.)/(K+l>)) 


LATERAL CEMTER, ,QF„ PRESSURE OF TH£ . R.I.QHT.J1ALF 

CPP • The spanwlse center of pressure of the 

right half of the wing nondlmenslonal Ized 
by BREF 

» CMBP/CLP 
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LATERAL , CENTER OF PRESSURE OP THg LgPT HAi P 

CPM - The spanwise center of pressure of the 
left half of the wing nondlmenslonat Ized 
by BREP 

• -CMBM/CLM 

PAJmiAa MOMENT 


CMP 


where 

CBAR 

and 

X 


• •*1*/(S*CBAR> 


If 

WING 


ACp*X*dS 


■ The reference chord/ 

• the dimensional chordwtse coordinate. 

* xsl^b/2 (see fig. 1) 


Although the above Integral can be evaluated 
analytically In the chordwtse direction/ It can 
not be analytically evaluated spanwise. The 
spanwise integration Is done using SUBROUTINE 
INTGRT. To allow the. user to assess the 

convergence this Integration Is done repeatedly 
with more and more points. The maximum number of 
points used Is the minimum of JJMAX and JJ. JJMAX 
Is the maximum number of available Integration 
points and Is read from the geometry file while JJ 
Is a number entered by the user. 

CHQRDWISE CEN TER OF PRESSURE 

XCP • The chordwtse center of pressure 
riondliTienslona 1 1 zed by CBAR 

- -CMP/CCL. 


10 



normal loads program 


4 USERlS INSIRUCtLQNS 

4.1 INITIAL SETUP FOR AMES‘ TSS SYSTEM 

For either batch or conversational processing the 
following TSS commands must be given. These commands 
are required once and only once for each user 10. The 
first three commands create the identification number 
file named iOFILE. This file contains four zeroes in 
binary form. 

SHARE MEDAN^FSARTM. INIDFILE 
CDS MEDAN, IOFILE 
DELETE MEDAN 

SHARE MEDAN, FSARTM,LSPROG. VI 

4.2 CONVERSATIONAL USE ON AMES* TSS SYSTEM 

All integer data should be entered in a 1615 
format, all floating point data in 8F10.0 
format^ and all logical data in lOLl format. 

USER: After logging on enter the following: 

AMES USYSLIB 
JOBLIBS SYSULIB 
JBLB MEDAN 

It is not necessary to issue DDEFs for anything 
except the input data since the program 

automatically issues them using the subroutines 
GEMFIL and BCFtL. 

USER: CALL FORCED 

PROG: ENTER BATCH 

USER: Enter carriage return for conversational mode. 

PROG: ENTER OOlSK (NEG. HALTS ) 

USER: For terminal output enter carriage return. For 

output to a disk file enter a positive non-zero 
number less than 10* Such a disk file will be 
referred to as the output file. For the AMES* 
TSS system the output will be found on the file 
named OUTPUT. FOR .Njt where X Is the numerical 
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value of ODtSK. The program Issues Its own DDEF 
commands so no control cards are needed. The 
value entered must be different from previous 
values for which the corresr>ondl ng output 
datasets have not yet been printed. The program 
uses logical unit 4 for this output. If a 
negative value Is entered/ the program will 
terminate. 

PROGt ENTER IDl/ ID2/ ID3/ ID4 

USER: Enter Identification numbers 

IDl Identification number of the geometry file 
from which the AIM and BC flies have been 
derl ved. 

ID2 Identification number of AIM file. 

1 03 Identification number of the boundary 

condition file. 

104 Identification number of the solution file. 

PROG: ENTER COMBINATION CODE 

USER: Enter combination code (LCOMB) to define a new 

set of weights as described below. If the 
output Is being placed on an output file (t.e. 
l^ODISKi^B)/ then the weight matrix will be 
printed on the output file. 

If LCOMB>0/ then LCOMB equals the number of 
combinations and the weights of each combination 
must be entered by the user. This data Is 
prompted for and entered In the following 
fashion: 

PROG: COMBINATION 1 

PROG: ENTER WEIGHTS OF FIRST NSYM CASES 

USER: Enter the set of weights. During 

execution the actual value of NSYM/ which 
Is the number of symmetric cases on the 
solution file/ is Inserted In the above 
message. After the above message Is 
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given the user should enter the NSYM 
weights corresponding to the symmetric 
cases of combination 1. 

PROG; ENTER WEIGHTS OF LAST NASYM CASES 

USER: During execution the actual value of 
NASYM, which Is the number of 

antisymmetric cases on the solution file. 
Is Inserted In the above message. After 
the above message Is given the user 
should enter the NASYM weights 
corresponding to the antisymmetric cases 
of combination 1. 

PROG: COMBINATION 2 

PROG; ENTER WEIGHTS OF FIRST NSYM CASES 

USER: Continue entering weights as above until 

the weights of all the combinations have 
been entered. 

If LCOMB«0, then the number of combinations 
equals the number of cases (solutions) and each 
solution with a factor of 1.0 Is treated as a 
combination, l.e., the weight matrix will be 

automatically computed as an Identity matrix. 

If LC0MB“-1, then the number of combinations 
equals the number of symmetric cases and each 
symmetric solution with a factor of 1.0 Is 

treated as a combination, l.e., the weight 
matrix Is an identity matrix with the last NASYM 
diagonal elements set to zero* 

If LC0MB“-2, then the number of combinations 
equals the number of antisymmetric cases and 
each anti symmetric solution with a factor of 1.0 
Is treated as a combination, l.e. the weight 
matrix Is an Identity matrix with the first NSYM 
diagonal elements set to zero. 

If LC0MBi.-3, then the current set of weights 

will be used. Do not enter -5 unless a set of 

weights has been previously entered. 
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PR062 * 

USER: At this point the user must begin entering 

commands. After each command he will be 
prompted for supplemental Input or with another 
* sign Indicating that he should enter another 
command. A full description of the commands Is 
given In section 5. After a NEW or KNEW command 
the next Input will be ODISK followed by IDI^ 
102/ 103/ and 104 followed (for the NEW command) 
by LCOMB followed by more commands. 

4.3 AMES’ TSS BATCH JOBS 

The batch mode operates the same as the conversational 
mode with the sole exception that a ”T“ must be put In 
column 1 on the first card. This “T'* suppresses all 
subsequent conversational prompts. 

4.4 OTHER COMPUTERS 

Remove all calls to 6EMFIL/ BNKFIL/ OBEY/ and CVRT In 
the main program and SUBROUTINE SBLOAD and use 
appropriate tape or disk control cards In. their place. 
These/ hopefully/ are the only changes that need to be 
made since considerable effort was made to code the 
program In standard FORTRAN. Then follow the 
instructions/ where appropriate/ In sections 4.2 and 
4,5. 
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5 DESCRIPTION QE fiOMMAHCiS 

A description of the commands^ which control the flow 
of the program^ is given in this section, in all cases 
the first three tetters of a command are sufficient 
input. Ait Integer data should be entered in a 1615 
format and alt floating point data in 8F10.D format. 
The inputs if any, associated with each command is to 
be entered on the fot towing tine In conversational 
processing and in batch processing is to be on cards 
immediately following the command card. Any input 
required is prompted, for in the conversational mode. 
The same command may occur more than once. This Is 
useful in correcting data entered in error. 

The commands which are the most basic and useful are: 
SPANLOADS, which causes the program to compute and 
print the span loading information; NETLOADS, which 
causes the program to compute and print the overall 
force and moment data; PRESSURE, which causes the 
program to compute and print pressure coefficient 
distributions; NEW or KNEW, which cause the program to 
start a new case; and STOP, which terminates execution. 
These and all the other commands are fully explained 
below: 


EFFECT: Causes the program to continue execution in 
the batch mode even if an invalid command is 
encountered. 

£££ 

INPUT: NTYPE plus other Input which varies with 

NTYPE. Valid values for NTYPE are 0, 1, 2, 3, 
and 5. 

EFFECT: Defines spanwise (ETA) stations at which the 

chordwise pressure distributions wilt be 
computed and printed following a PRESSURE 
command. If this command Is not given, the 
spanwise control points on the geometry file 
will be used. This command allows a direct 
comparison to be made with other theories 
and/or experimental data. The input required 
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for the various values of NTYPE and the 
spanwlse stations defined thereby are given 
below: 

INPUT FOR.NTYPE-0: NONE 

The following type of stations are 
prepared: 

ETA-COS(i*Pi/(MREF+l)) for l.l^MM. 

MREF comes from the solution file. The 
equation solving program transferred this 
number from the Influence matrix file. 
The influence matrix program obtains this 
number either from the geometry file or 
as user input after the MREF command. If 
obtained from the geometry file/ MREF 
will be equal to the geometry program 
variable NN, 

INPUT FOR NTYPE-lt NQ 

The following type of stations are 
prepared: 

ETA-C0S<I*PI/(JJMAX^1)) for l-MQ, JJMAX^NQ 

In the conversational mode the program 
prompts the user for the value of NQ, 
The maximum allowable value for NQ Is 
equal to the value of NDfM3 in SUBROUTINE 
SBLOAD. Currently NDIM3-47. JJMAX is 
the number of spanwlse stations at which 
the wing data is given. JJMAX ts read 
from the geometry file. 

INPUT FOR NTYPE-2: NSTA 

The following type of stations are 
prepared: 

ETA-COS(i*PI/(NSTA+l)) for i-I,NSTA. 

in the conversational mode the program 
prompts the user for the value of NSTA. 
The maximum allowable value for NSTA is 
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equal to the value of N0IM3 tn SUBROUTINE 
SBLOAD. Currently NOIM3-47. 

INPUT FOR NTYPE-5: a table of spanwise (ETA) 

stations 

In the conversational mode the program 
prompts the user to enter the table. 
This table Is to be entered one value per 
line or card. Values entered need not be 
In any particular order. The last entry 
must be followed by a line or card 
containing a number greater than 1.0 In 
order to signal the end of the table to 
the program. The allowable number of 
entries Including the value greater than 
1*0 is equal to the value of NDIM3 tn 
SUBROUTINE SBLOAD. Currently NDIM3-47. 

INPUT FOR NTYPE-5J ETMI N, ETMAX^DETA 

The following type of stations are 
prepared : 

ETA»ETMIN/ETMIN+0ETA^ETMIN+2*DETA^ ... ^ 
ETMAX. 

ETMAX and DETA have default values. The 
default for DETA Is ETMIN^ while the 
default for ETMAX Is 1.0. The user must 
Insure that no more than NDIM3 stations 
are defined In this way. N0IM3 Is a 

variable tn SUBROUTINE SBLOAD and 
currently Is equal to 47. 


ETAS 

INPUT: Same as for ECP command 

EFFECT: Defines spanwise (ETA) stations at which the 

various spanwise loading distributions will be 
computed and printed following a SPANLOADS 
command. If this command is not glven^ the 
spanwise control points on the geometry file 
will be used. This command allows a direct 
comparison to be made with other theories 
and/or experimental data. Follow the 
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Instructions for the ECP command with the 
exception that the maximum number of stations 
Is determined by NDIMl In SUBROUTINE SBLOAD. 
Currently NDIMl-200. 


KNEW 

EFFECT: Performs the same function as NEW (see below) 

with the exception that the current stations 
for pressures and spanloads are retained. 


HETIQAD^ 

INPUT: JJ.NOUT 

EFFECT: Computes' and prints the following overall 

results: CCL (lift coeff 1 dent) # CMP (pitch 

moment)^ CMR (roll moment)# left and right 
root bending moments# left and right lifts# 
left and right centers of pressure# CDI# and 
the vortex drag factor. JJ Is the maximum 
number of points to use In integrating for the 
pitching moment. The default for JJ Is JJMAX# 
which comes from the geometry file. If JJ 
exceeds JUMAX the program will use only up to 
JJMAX points. NOUT Is the output level for 
subroutine INTGRT (see listing). NOUT»0 Is 
the usual choice, 

NEW 

EFFECT: Starts a new case. After this command Is 

given the program returns to the point at 
which ODISK Is requested (section 4,2). The 
value entered for ODISK must be different from 
previous values for which the corresponding 
output datasets have not been printed. 
Subsequently IDl# ID2# ID5# and ID4 are 
requested. Then# after the user enters the 
new set of IdentlHcatlon numbers# the program 
reads a new solution file and# If necessary# a 
new geometry file and resets the spanwise 
stations at which the spanwise loading 
distr Ibutlons are calculated and the spanwise 
and chordwise stations at which the lifting 
pressure distribution Is calculated. Then the 
program requests the combination code (LC0M6) 
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and/ If LCOMB>0/ a new set of weights as 
described in section 4.2. Following this the 
program is In the command mode again. 


EFFECT: Computes and prints the lifting pressures at 

the chordwise locations defined by the XCP 
command and at the spanwise locations defined 
by the ECP command* if the XCP or ECP 
commands have not been given/ then the 
chordwise and spanwise stations used will be 
the spanwise control points. 

PWEL6H.T.S 

INPUT: NONE 

EFFECT: Prints the current weights matrix on the 

user's terminal. This command is generally 
only used In the conversational mode. 


&PAHLQADS 

EFFECT: Prints the loading distributions at the 

spanwise stations defined by the ETAS command, 
if the ETAS command has not been given/ then 
the spanwise control points will be used. 

HOC 

EFFECT: Halts execution. 

TSS 

INPUT: A TSS command of 80 characters or less. 

EFFECT: The command Is passed to the AMES' TSS 

operating system. After the system processes 
the command/ control returns to the program. 
This command Is a special one for the AMES' 
TSS version of the program, 

WEIGHTS 

INPUT: LCOMB and/ if LCOMB>0/ a new set of weights. 
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EFFECTt Allows the user to define a new set of weights 
In the same manner as when starting a new case 
(see section 4*2). 

XCP 

INPUT} NTYPE plus other Input which varies with NTYPE 

EFFECT} NTYPE defines the chordwise (x) stations at 
which the lifting pressures will be computed 
following a PRESSURE command. The quantity x 
Is the local chordwise variable such that 
-l^xi+l. If this command Is not glven^ the 
stations will be Identical to the spanwise 
control points on the geometry file (using the 
correspondence x«eta}. This command allows a 
direct comparison to be made with other 
theories and/or experimental data. The Input 
required for the various values of NTYPE and 
the chordwise stations defined thereby are 
given below} 

INPUT FOR NTYPE-0} NONE 

The following type of stations are 
prepared} 

x-COS(l*Pl/(MREF + D) for 1-1, MM. 

MREF comes from the solution file. The 
equation solving program transferred this 
number from the Influence matrix file. 
The Influence matrix program obtains this 
number either from the geometry file or 
as user Input after the MREF command. If 
obtained from the geometry file, MREF 
will be equal to the geometry program 
variable NN. 

INPUT FOR NTYPE-1: NQ 

The following type of stations are 
prepared} 

x-COS(l*PI/(JJMAX+D) for 1 -NQ,dJMAX,NQ. 
In the conversational mode the program 
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prompts the user for the value of NQ» 
The maximum allowable value for NO Is 
equal to the value of NDIM2 In SUBROUTINE 
SBLOAD. Currently NDtM2«100. JJMAX Is 
the number of spanwise stations at which 
the wing data is given. JJMAX ts read 
from the geometry file. 

INPUT FOR NTYPE-2: NSTA 

The following type of stations are 
prepared: 

x-COS(l*PI/(NSTA+D) for l-l^NSTA. 

In the conversational mode the program 
prompts the user for the value of NSTA. 
The maximum allowable value for NSTA ts 
equal to the value of NO I M2 In SUBROUTINE 
SBLOAD. Currently NDIM2«100. 

INPUT FOR NTYPE-3: table of values of x 

In the conversational mode the program 
prompts the user to enter the table. 
This table ts to be entered one value per 
line or card* Values entered need not be 
in any particular order. The last entry 
must be followed by a line or card 
containing a number greater than 1.0 in 
order to signal the end of the table to 
the program. The allowable number of 
entries Including the value greater than 
1.0 is equal to the value of NDIM2 in 
SUBROUTINE SBLOAD. Currently NDIM2>100. 

INPUT FOR NTYPE"4: table of values of CHI 

(CHI-(x+l)/2.) 

In the conversational mode the program 
prompts the user to enter the table. 
This table Is to be entered one value per 
line or card* Values entered need not be 
In any particular order. The last entry 
must be followed by a line or card 
containing a number greater than 1.0 in 
order to signal the end of the table to 
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the program* The allowable number of 
entries Including the value greater than 
1.0 is equal to the value of ND1M2 In 
SUBROUTINE SBLOAD. Currently NOIM2*100. 

INPUT FOR NTYPE-5: XMIN^XMAX.DX 

The following type of stations are 
prepared: 

x»XMIN,XMIN>DX^XMIN<^2*DX, ... , XMAX. 

XMAX and OX have default values. The 
default for XMAX Is 1.0/ while the 
default for DX Is XMIN-C-I.). The user 
must Insure that no more than NDIM2 
stations are defined In this way. NDIM2 
Is a variable In SUBROUTINE SBLOAD and 
currently Is equal to 100. 
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6 SAMPLE TERMINAL &£aS.,LQIi 

A sample conversational terminal session on the Ames' 
360/67 TSS computer system Is reproduced In this 
section with additional comments added In parentheses. 
During this session the 5 sets of solutions obtained In 
the sample case of the equation solving program (ref. 
U) were utilized. These sets of solutions were for an 
aspect ratio 2 , rectangular wing and differed from one 
another only In the number of control points and/or 
pressure modes used. 

The first set (IDit«15) was first used to Illustrate the 
simplest treatment^ l.e»# only the NETLOADS and 

SPANLOAOS commands. The default weight matrix and 

default spanwise stations (the spanwlse control points) 
were used. Next the use of the WEIGHTS command is 
Illustrated. Then the chordwise and spanwlse stations 
for the lifting pressure distribution were set up and 
the pressure distribution obtained. These stations 
were chosen to be Identical to the stations used In 
ref. 6 so that a direct comparison among the 4 theories 
could be made. Then the spanwlse stations used In ref, 
6 for the spanwlse load distributions were set up (NSTA 
■ 15 ■ the number of spanwlse control points used In 

ref. 6) and the spanwlse loading was obtained. This Is 
all that was done with the first solution set. 

For the remaining 4 solution sets the various stations 
existing at the end of the first solution set and the 
weights matrix were retained (using the KNEW command 
and LC0MB"-3) In order to facilitate a direct 
comparison among the sets of results obtained using the 
author's computer programs and those obtained In ref, 
6. For each of these sets net loads, span load 
distributions, and pressures were obtained. 

The output from the following test run was directed to 
disk files, which were later printed and are given In 
appendix I : 

LOGON user Id, password, termi nal Id 

AMES USYSLIB 

JOBLIBS SYSULIB 

JBLB MEDAN 

DDNAME-JBLBOOOl 

CALL FORCES 
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ENTER BATCH 
(carriage return) 

ENTER ODISK (NEG. HALTS } 

1 

OUTPUT IS ON OUTPUT. FOR. N1 
CANCELLED: DDNAME FT04F001 UNKNOWN 

(Messages such as the above occur because of the 
automatic file defining feature of the Ames' version. 
They do not indicate any error.) 

ENTER IDl, ID2, ID3, ID4 
4 2 9 15 

CANCELLED: DDNAME FT07F001 UNKNOWN 
CANCELLED: DDNAME FT12F001 UNKNOWN 
ENTER COMBINATION CODE 
0 

♦ 

TSS 

CPUTIME? 

3.796 SECONDS 

♦ 

NETLOAOS 

ENTER JJ AND NOUT 
0 0 

♦ 

SPAN LOADS 
+ 

WE 1 6HTS 

ENTER COMBINATION CODE 
1 

COMBINATION 1 

ENTER WEIGHTS OF FIRST 2 CASES 

1 . 1 . 

ENTER WEIGHTS OF LAST 1 CASES 

1 . 

♦ 

SPANLOADS 

♦ 

NETLOADS 

ENTER UJ AND NOUT 
50 

* 

ECP 

ENTER NTYPE 
3 

ENTER TABLE OF ETAS— 1 PER LINE ENDING WITH VALUE 
GREATER THAN 1 
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0 , 

0.3827 
0. 7071 
0.9239 
1 
+ 

XCP 

ENTER NTYPE 
4 

ENTER TABLE OF CHI VALUES 

0.005 

0.0125 

0.025 

0.05 

0.1 

0.15 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

0. 95 
1 

♦ 

WEIGHTS 

ENTER COMBINATION CODE 
1 

COMBINATION 1 

ENTER WEIGHTS OF FIRST 2 CASES 

1 . 

ENTER WEIGHTS OF LAST 1 CASES 

0 . 

♦ 

PRESSURES 

ETAS 

ENTER NTYPE 
2 

ENTER NSTA 
15 

+ 

SPANLOADS 

+ 
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KNEW 

ENTER ODISK (NEG. HALTS > 

2 

OUTPUT IS ON OUTPUT. FOR. N2 
ENTER IDl, ID2, ID3, ID4 
4 2 9 16 

ENTER COMBINATION CODE 

-3 

♦ 

NETLOADS 

ENTER JJ AND NOUT 
50 

♦ 

SPANLOADS 

♦ 

PRESSURES 

+ 

KNEW 

ENTER ODISK (NEG. HALTS ) 

3 

OUTPUT IS ON OUTPUT. FOR. N3 
ENTER IDl, ID2^ IDS, ID4 

4 2 9 17 

ENTER COMBINATION CODE 

-3 

+ 

NETLOADS 

ENTER JJ AND NOUT 
50 

+ 

SPAN LOADS 
+ 

PRESSURES 

+ 

KNEW 

ENTER ODISK (NEG. HALTS ) 

4 

OUTPUT IS ON OUTPUT. FOR. N4 
ENTER IDl, ID2, ID3, ID4 
4 2 9 18 

ENTER COMBINATION CODE 

-3 

+ 

NETLOADS 

ENTER JJ AND NOUT 
50 
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SPANLOADS 

+ 

PRESSURES 

+ 

KNEW 

ENTER ODISK (NEG, HALTS ) 

5 

OUTPUT IS ON OUTPUT. FOR. N5 
ENTER IDl, ID2, ID3, ID4 
4 2 9 19 

ENTER COMBINATION CODE 
-3 

+ 

NETLOAOS 

ENTER JJ AND NOUT 
50 

+ 

SPANLOADS 

+ 

PRESSURES 

+ 

TSS 

CPUTIME? 

17.438 SECONDS 

+ 

STOP 

TERMINATED: STOP 

PRINT OUTPUT.FOR.Nl^PRTSP-EDIT^STATION-RMTOS 
PRINT BSN-9132^ 200 LINES 

PR I NT OUTPUT. FOR. N2, PRTSP»ED I T^ STATI ON-RMTOS 
PRINT BSN»9133^ 100 LINES 

PRINT OUTPUT. FOR. N3^PRTSP-EDIT^STATION-RMT05 
PRINT BSN«9134^ 100 LINES 

PRINT OUTPUT. FOR . N4^ PRTSP-ED 1 T, STAT I ON-RMTOS 
PRINT BSN-9135, 100 LINES 

PRINT OUTPUT. FOR. N5^ PRTSP-EDIT#STAT!ON«RMT05 
PRINT BSN-9136^ 100 LINES 

LOGOFF 
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LOADS ON THIN, LIFTING WINf, 

SSSBCCBSSSSmBBSCBSBSSCSSSSs 


rectangular wing AR b 2 


IDI 

e 

4 

10? 

8 

? 

103 

ii» 

m 

R 

lOtt 

8 

IS 

OELTAO 

e 

«,0000 

EPS 

m 

0,5000 

MACH 

8 

0,0000 

JJ 

8 

l9l 

jjhax 

a 

IRI 

NMAX2 

9 

5 

KH2 

a 

11 

PPNfW 

a 

5 

MHNEH 

8 

n 

NROWS A 

a 

3 n 

CWTYPE 

8 

0 

SWTYPC 

a 

0 

RCS 

9 

TTFFFFFFFF 

BCAS 

s 

FTFf FFFFFF 

hACm 

(B/2)/BPEF 

B O’. 0 0 00 

s 1,0000 

CSAfl/BREE 

a 1,0000 

ASPECT 

RATIO s Z'.OOOO 


J 1 p i 3 • 7 ^ 


OUTPUT. FOn.Nl 



weights 


COHBINaTION j 

r.OOOOOO 0',000000 o'.oooooo 
COMBINATION 2 

0,0000 00 r.OOOOOO o'.oooooo 
combination s 

o'.oooooo o'.oooooo r.OOOOOO 
combination j 

CL ■ 2,470170 

induced drag. 9 0.970900 

VORTEX DRAG FACTOR ■ 1,000606 

ROLLING MOMENT M/ (0*S*2*BKEF) m 0,000000 

Right root bending moment MBR/(QtiS«BREFl • 0, 529675 

left root bending moment MBL/(0*S*BRFf ) ■ 0,529675 

LIFT ON ETA.GT'O b 1.237068 

LIFT ON eTA.LT-,-a - -M 1.237066 

CENTER OF PRESSURE OF RIGHT HALF, V/BRFF s O’,O20163 

center of pressure of Left half, y/brff • -o,o2ei63 

pitching moment about X b 0 and center of pressure 

J CM/(G*S*CBaR) X(C,P,)/CBaR 

11 •0,518157 0,209026 

JiS ;0^161&7 Q.2i)9oeAi 

07 1,0,516157 0.209026 

95 •0, 516157 0.209026 

191 »0, 516157 0.209426 



COMBlNiTION ? 


CL s l,<?b6263 

induced drag a 0,610638 

VORTEX drag factor a 1,00355? 

POLUII^G MOM6NT M/{0*8*|ApREEi a 0,000000 

RJGHT ROOT PENDING MOMENT HRfi/ (Q*S*BRFF 1 a 0,«22256 

left root pending moment mRL/(T*S*BREF) a 0,422256 

LIFT ON ETA.GT 0 a 0,978132 

lift on ETA.LT.O b 0,978132 

center of pressure OF RIGHT HALF, V/RrEF = 0‘,«5J69T 

center of pressure of Left half, v/bref s -o.usisr? 

pitching moment about X a ft AND CENTER OF PRESSURE 


J 

CM/(0*S*C84R) 

X(C.P, )/C0AR 

It 

-0,T56382 

0,386646 

25 

-0,756582 

0,386646 

4? 

«0. 756562 

0 386646 

95 

•0,756562 

0,586646 

191 

-0,756582 

0.586646 

combination 

3 

--- 



CL ■ 0,000000 

induced drag 3 0,183247 

rolling moment H/(0*S*2*RREF ) a 0.189680 

RIGHT ROOT bending MOMENT mrR/ ( 0«S*rREF ) a 0,189680 

left root bending moment MBL/ {D*S*BREF> a -0,189600 

LIFT ON ETA.GT.O a 0,32Ul9 

LIFT ON ETA.LT.O B -0, 521319 

CENTER OF PRESSURE OF RIGHT HALF, Y/BPEF a 0,590318 

center of pressure of left half, Y/BREP a -0,590318 

pitching moment ABOUT X a ft AnD CEnTEr OF PRESSURE 
J CM/(0*S*CBaR) X(C,p,)/C8aR 

11 o'oooooo 

25 0,000000 

47 0,000000 

95 0,000000 

191 0.000000 



FTa 

CL#C/2s 

CL 

CL*C/CCl*CAVG 

2CMfi/4)C/W C,P,fl/«) 

COHBINaTIOn 

1 




0,^65926 

0,866025 

0,707107 

0,500000 

0,258819 

0,000000 

0;211735 

0,«05291 

0,565220 

0,682165 

0;752i»95 

0.7758«a 

n',fl969«l 

1,621162 

2,260680 

2,728661 

3,009970 

3,103375 

o' 342312 
0,655233 
0,913791 
1,102856 
1,216554 
lU^9306 

0^067269 pO, 079426 
0,105695 •0,065197 
0,113871 •0,050366 
0,106477 -0,039022 
0,097121 -0,032266 
0,093296 •0,030065 

C0M9INATI0N 

? 




0;96S926 

0,866025 

0;707107 

0,500000 

0,258819 

0,000000 

o' 175510 
0,329908 
0‘-«5207« 
0;S38337 
0,509050 
0.605716 

0. 702039 
1.319633 

1 , e0829fl 

2, !535«8 
2,556199 

2,a2286a 

o' $58867 
0,674560 
0,924361 
1,1 00744 
1,204437 
1^38515 

•0,109422 0.155865 
•0,196098 0,148601 
•0,255870 0,141496 
«0,29265« 0,155907 
•0,311984 0,132410 
•0,317951 0,131230 

combination 

5 




0,965926 

0,866025 

0,707107 

0,500000 

0,258819 

0,000000 

0' 121606 
0,210255 
0,2616a9 
0,208073 
0,119579 
0.000000 

0,686431 

0,841018 

0,966596 

0,832293 

0,478316 

0,000000 

•0,000000 

•0,000000 

•0,000000 

•0,000000 

•0,000000 

•0.000000 

0,066301 -0,136301 
0,099890 •0,118773 
0,096968 •0,100319 
0,07l«47 .0,085844 
0,036796 -0,076927 
0,000000 .0,077357 


WEIGHTS 


COMBIHaTION 1 

1,000000 l ',000000 l ’,000000 



e t4 

CL*C/2b 

CL CL*C/CCL*CAVG 

2CMft/9)C/0 

C,P, f 1/91 

combination 

\ 





0,9659?6 

0' 508855 

2,055912 

0^956915 

0',024199 

.0,011669 

C,B(>bO?5 

0,905«5S 

3.701612 

0,85.5597 

0,009906 

•0,002506 

O,70TJ07 

|'256936 

5,055799 

1,136625 

-0,045031 

0,008992 

0,500000 

0,258819 

0,000000 

i:«2856fl 

5,719272 

1,209779 

-0,119729 

0,020078 

1,961112 

5,699998 

1,319157 

•0,176067 

0,030968 

1,381555 

5,526218 

1,297328 

-0,229656 

0,090653 

-0,258819 

1,221951 

9,807802 

i;i03230 

•0,251658 

0,051987 

«o,?ooooo 

1,012922 

9,049686 

0,919059 

.0,257629 

0,063616 

•0,707107 

0,775695 

3 102578 

0,700287 

-0,230967 

0,077022 

•0,866025 

0,529999 

2.099776 

0,973993 

•0,190293 

0,090626 

•0,965926 

0.265637 

1,062597 

0U3<>829 

.0,108959 

0,102069 

combination 

a M A A 





CL 

^ w * • 

8 

9,930935 




induced drag • 

3,311915 




VORTEX DRAG 

FACTOR • 

1.059986 




rolling moment m/(Q*8*?*R«FE1 

s 

0,1 89660 


RIGHT ROOT 

bending moment m8R/(0*S*BrFF) b 

1.191611 


left ROOT 

BENDING MOMENT H«L/ f O aS*0REF ) 8 

0,762251 


LIFT ON ETA 

,GT.0 b 

2,556537 




LIFT ON ETA 

.LT.O * 

1,093900 




CENTER OF PRESSURE OF RIGHT HALF 

, T/BPEF s 

0'.950067 


center of pressure of 

LpFT half 

, Y/BReF 8 

-0,902977 


PITCHING MOMENT ABOUT X 

a 0 and 

CENTER OF PRESSURE 


J 

CM/(0*8*CBaR) 

X(C,p.)/CBaR 



1 1 

-1,279555 


0,287677 



23 

-1,279535 


0,287677 

0,267677 



A7 

-1,279555 






WElfiNTS 


combination I 

liOOOOOO o'oooooo o'oooooo 
combination 1 


oloooooo 


ETA s 

•■ 0,990000 
• 0,975000 
• 0,950000 
• 0,900000 
• 0,000000 
• 0,700000 
• 0,600000 
«o,aooooo 
• 0,200000 
• 0,000000 
0 200000 
0,400000 
0^600000 
O.AOOOOO 
0,900000 


0* 005000 
0,012500 
0,025000 

o;o 9 oooo 

0,100000 

0,150000 

0,200000 

0,300000 

0,400000 

0,500000 

0,600000 

0,700000 

0,000000 

01900000 

0^50000 


OELTa*CP 
311516403 
191775360 
13,798651 
9,498217 
6,356733 
4;904237 
4^006256 
2,094391 
2,198984 
1,706087 
I , 328721 
1,020911 

01751140 

0;485060 

01329260 



0'.3flJ>700 


ETA « 


X 

CHI 

delta. CP 
so' 069539 

•0,990000 

o' 005000 

•0,975000 

0,012500 

18,851593 

•0,950000 

0' 025000 

I3;i351t8 

•0,900000 

0,050000 

9 019387 

•0,800000 

oj lOOOOO 

5,999567 

«0,700000 

0,150000 

9,599013 

•0,600000 

0,200000 

3,727965 

■o,aooooo 

ojiooooo 

2,658195 

•0,^00000 

0' <100000 

1,997970 

•0.000000 

o,aooeoo 

0,500000 

1;537669 

0,600000 

1,192592 

0.000000 

0,700000 

0,915512 

0 600000 

0 600000 

0, 673715 

0,800000 

0J900000 

0.939571 

0,900000 

0,950000 

0.2992«9 

ETA * o’. 

X 

707100 

CHf 

oelta.cp 

*0,990000 

0,005000 

25.631666 

•0,975000 

0,012500 

15,983666 

•0,950000 

0,025000 

11,039855 

7,999990 

*0,900000 

oCosoooo 

•0,800000 

0 ; lOOOOO 

9^800809 

•0,700000 

0,150000 

5,577179 

•0,600000 

0J200000 

2,831 959 

•0,^00000 

0 300000 

1,991271 

•0,200000 

0,900000 

1,921001 

•0,000000 

0,500000 

1:076«85 

0,200000 

0,600000 

01832978 

0,U00000 

0,700000 

0,6599?6 

0,600000 

0, 600000 

0,972996 

0,600000 

0,900000 

0,306167 

0,900000 

0,950000 

0,207928 



ETA B 0’,9S3900 


K 

CHI 

OtLTAfCP 



•>0,990000 

0* 005000 

171776794 



•0,975000 

0,012500 

10,747533 



•0,950000 

0,025000 

7,059169 



•0,900000 

0,050000 

4,333611 



•0,800000 

0,100000 

2,389219 



•0,700000 

0, 150000 

1,614014 



•0,600000 

0,200000 

1,239967 



•0,400000 

0,300000 

0,912861 



•0,200000 

0,400000 

0,759190 



•0,000000 

0,500000 

0,611839 



o;2ooooo 

0,600000 

0,445402 



0,400000 

0,700000 

0,286298 



0,600000 

0,800000 

0,174114 



0,600000 

0,900000 

0,127857 



0,900000 

0;950000 

0,113100 



eta 

CL*C/2b 

CL CL*C/rCL*CAVS 

2CMf 1/41C/8 

C 

combination 

! 




0,900785 

0' 159810 

0,639292 0'?S836S 

0,052435 

• 0 

0,923860 

0,311866 

1,247465 0,504194 

0,090029 

• 0 

0,851470 

0,448907 

1,795628 0,725747 

0,110096 

• 0 

0,707107 

0,565220 

2,260880 0,913791 

0, U3871 

• 0 

0,555570 

0,657235 

2.620938 1,062550 

0,108945 

• 0 

0,562681 

0;721222 

2,892808 1,169231 

0,101427 

• 0 

0 195090 

O' 762717 

3,050869 r?530fl4 

0,095487 

-0 

0,000000 

0.775844 

3,103375 1'254306 

0,093296 

• 0 


P,M/4j 


082027 

072SU 

061313 

OS0366 

041441 

035061 

031298 

030063 



LOADS OM thin, lifting WINfi 
BB>«sssssssszsstessaesssssssc 


PECTANGULA 9 t»JlNG AH a 3 


101 

s 

4 

ID? 

n 

? 

IDS 

8 

P 

I0« 

8 


DELTAO 

8 

4,0000 

EPS 

8 

0,5000 

HACH 

8 

0.0000 

JJ 

8 

I’l 

JJHAX 

O 

IPI 

NHAX2 

8 

5 

KK2 

8 

5 

PPNEta 

8 

5 

MHNeW 

8 

11 

NPOWSA 

8 

So 

chtype 

8 

0 

SWTYPE 

i 

0 

BCS 

8 

TTFFFf FFFF 

BCAS 

D 

FTFFFf FFFF 

HACH 

(e/2)/0REF 

» 0>000 
■ 1,0000 
• I.OOOO 

CSAR/BREF 

ASPECT 

PATIO a 2.0000 


ntiTPUT.rnR 



WFIGHTS 


■p- 

3 


COW81NATION I 

r.OOOOOO O'.OOOOOO 0* 000000 


combination 


CL ■ S,aTl5tt9 

INDUCED DRAG ■ 0,972838 

VOPTEX DRAG FACTOR a l, 000 fe<l 8 

rolling moment ^ p 

R|OmT root BENOiNG moment M 0 R/ (Q*S*HREF) a 0 *S 

left root bending moment Mpi./rQ.S* 0 R£F » - 0*5 

lift on ETA.GT.O a 1,235778 * 

LIFT ON ETA.UT.O a J , 235776 

CENTER OF PRESSURE OF RIGHT HALF, V/BREF a 0 * «, 

pressure of LfFT half, V/Bwpf b *o!a( 

PITCHING MOMENT ABOUT X a n AND CENTER OF PRESSURE 


CM/{(Jts*CBAR) 
•0,518607 
aO, 518607 
.0,516607 


X(C,P,1/C8aR 

0,209750 

0,209750 

0.209750 


eta 

CL*C/2R 

CL 

cl.c/ccl*cavg 

2CMf 1 /6)C/P 

C,P,f 1/6) 

combination 

! 





0,980785 
0,923880 
0,831670 
0,707107 
0,555570 
0,382683 
0 195090 
0,000000 

■ ^ 

0* 159771 
0*311398 
0' 668387 
0:566585 
0,656567 
O’ 722670 
0,761899 
0:776995 

0,639085 

1,266393 

1,795569 

2,258339 

2,626186 

2,869881 

3,067595 

3,099979 

0* 256576 
0,506296 
0,725677 
0,913733 
1,062566 
1,169257 
i;23SO60 
i'256266 

0,052376 
0,09l2bO 
0,110266 
0,112395 
0,106105 
0 099297 
0,095667 
0,096376 

•0,081956 

•0,073219 

•0,061689 

-0,069769 

•0,060603 

•0,036360 

-0,031319 

•0.030663 



combination 


1 


ETA a OlOOOOOO 


X 

CNf 

oei.ta#cp 

•0,990COO 

0‘, 005000 

31,872513 

•O.^TSOOO 

0,012500 

19,950165 

•0,950000 

0,025000 

13,668690 

•0,900000 

0 osoooo 

9,987066 

•0,B00000 

OMOOOOO 

6,297503 

•0,700000 

0,150000 

9,892097 

• 0, f^OOOOO 

0^200000 

3,956260 

•o,aooeoo 

0,300000 

2,877397 

•0,300000 

0,900000 

2,209339 

•0,000000 

0^500000 

1,717292 

0,200000 

0,600000 

1,333607 

0,000000 

0,700000 

1,016199 

OjbCOOOO 

0,600000 

0,792330 

oJbooooo 

0J900000 

0,962797 

0,900000 

0^950000 

0.332925 

ETA a o'. 

362T00 


X 

CHl 

DELTa.CP 

•0,990000 

o' 005000 

29’ 532990 

•0,975000 

0^012500 

18,578003 

•0,950000 

0,025000 

13,011273 

*0,900000 

0,050000 

9,006115 

•0,fl00000 

0, lOOOOO 

6,055693 

•0,700000 

Q , 150000 

9,661661 

•0,600000 

0^200000 

3,781092 

•0,000000 

o;3ooooo 

2,670871 

•0,200000 

0,900000 

1,9825u3 

.0,000000 

0^500000 

1,516655 

0,200000 

0,600000 

1,182029 

0,900000 

o;rooooo 

0,919961 

0,600000 

0,600000 

01686526 

0,600000 

0,900000 

0,991160 

0,900000 

0,950000 

0:295637 



FT* c 0, 

707100 


% 

CHI 

DIUTA^Cp 

•0,990000 

0' 005000 

25* 732437 

•0,975000 

0' 012500 

16;017560 

•0,950000 

0,025000 

0,050000 

11,032462 

•0,900000 

7,410202 

•0,000000 

OMOOOOO 

4,746520 

•0,700000 

0,150000 

5,533155 

■0,600000 

0,200000 

21802224 

•0,400000 

O' JOOOOO 

1,941195 

•0,200000 

0,400000 

1,437692 

•0,000000 

0,500000 

1,097511 

0,200000 

0,600000 

0,843628 

0,400000 

0,700000 

0,640770 

0,600000 

0,600000 

0,468129 

0,000000 

0,900000 

0,304315 

0,900000 

0'.950000 

0,209267 

ETA ■ 0, 

925900 


% 

rHi 

OEUTAiCP 

■0,990000 

0' 005000 

17',689465 

•0,975000 

0,012500 

0;029000 

10,809661 

•0,950000 

7,092859 

■0,900000 

0,050000 

4,3«aS?6 

«O,B0O0O0 

0,100000 

2,381964 

•0,T00000 

0,150000 

1,599733 

•0,600000 

0,200000 

1,216279 

■0,400000 

o; 3 ooooo 

0,898071 

•0.200000 

0,400000 

0,750157 

•0,000000 

0,500000 

0j609l06 

0,200000 

0,600000 

0J447497 

0,400000 

0,700000 

0,290559 

0,600000 

0,600000 

0,177410 

0,800000 

0,900000 

0,127898 

0,900000 

0^950000 

0;iU598 



LOADS 

ON THIN, lifting 

ssssss 

It 

It 

u 

li 

ssssssssasszs 

RECTANGULAR NIMG AR s 

IDI 

e 

H 

ID2 

m 

? 

103 

e 

9 

!D4 

1 

17 

DELTAO 

II 

ti.OOOO 

EPS 

a 

o.sooo 

MACH 

a 

0,0000 

JJ 

a 

1«1 

jjmax 

s 

191 

NMAX2 

a 

5 

KKi 

a 

% 

PRMEW 

a 

? 

MHNCW 

a 

0 

mROwSA 

a 

15 

cwtype 

a 

0 

shtype 

a 

ft 

0CS 

a 

TTf FFFFEFF 

nCAS 

a 

FTFFFFFFFF 

MACH 

(B/2J/BRCF 

• o'oooo 
■ 1,0000 

* 1.0000 

CSAR/BREP 

ASPECT 

Ratio s ?.oooo 


U*ll»73 


OUTPUT. FOR. M3 



weights 


combination j 

I ’,000000 o', 000000 o' 000000 
combination 1 


CL 

m 


,474199 




INDUCED DRAG a 

Oi 

,974926 




VORTEX DRAG 

FACTOR ■ 

1 . 

000648 




ROLLINS MOMENT m/(Q*s*?*rREF) 

9 

oloooooo 


riOhT root 

BENDING moment 

MBR/(Q*SaBREF) b 

0,529666 


lirt root 

bending moment 

M6L/(0*S*BREF) B 

0,529686 


LIPT ON ETA 

,0T^0 a 

1. 

257101 




LIFT ON ETA 

.UT « . 

w 

1. 

217101 




center of pressure or right half 

1 Y/BREE a 

01428166 


CENTER or pressure Or 

LEFT HALF 

, V/BREF b 

•0,426166 


PITCHING MOMENT ABOUT X 

m 

0 AND 

CENTER OF PRESSURE 


J 

CM/(Q*S*CBaR) 


x(C,p,)/C0aR 



11 

•o,sieteo 



0,209434 



2S 

•o.Biaieo 



0,209434 



til 

-0,B18180 



0,209434 



eta 

CL*C/2r 


CL 1 

CL*C/CCL*CAVr, 

2CMf 1/4)C/B 

C,P,fl/4) 

combination 

1 






0.9B0T85 

0'159VT2 

0, 

639886 

01256623 

0l050629 

.0,079122 

0,923880 

01311967 

1. 

247869 

01504352 

0,069493 

•0,071716 

0;S3U70 

0' 948661 

1. 

795525 

Oir25698 

0.110412 

-0,061493 

oJtotjot 

o'b6Sl74 


260697 

01913T07 

0,114814 

•0,050787 

0,55B570 

0;<»57223 


628692 

11062521 

0,109328 

•0,041567 

0-382S63 

01723231 


892922 

11169234 

0,101255 

-0,035001 

0,195690 

01762726 

3. 

050904 

11253086 

0,095298 

•0,031236 

0,000000 

01775848 

5. 

103393 

11254300 

0,093217 

-0,030037 



COHbTNaTION' 


1 


ETA • O'.OOOOOO 


K 

CNI 

delta. CP 

•0,990000 

0^005000 

5r,509599 

•0,975000 

0,012500 

19,77107? 

•0,950000 

0,025000 

13,795799 

•0,900000 

0,050000 

9,996516 

•0,400000 

0,100000 

6,356160 

•0,700000 

0,150000 

9,909528 

•0,h00000 

0,200000 

9,006797 

•o,(jooooo 

0' 300000 

2,895579 

•0,200000 

0,900000 

2,199999 

•0,000000 

0,500000 

1,706737 

0,200000 

0,600000 

1,3289<|9 

0,900000 

0,700000 

1,020793 

0,600000 

0,800000 

0,750909 

0,400000 

0,900000 

0,989978 

0,900000 0.950000 

TA a 0*,|a2700 

0^29289 

X 

CHI 

delta. CP 
30’,023071 

•0,990000 

O" 005000 

•0,975000 

0,012500 

18,827209 

•0,950000 

o;025000 

15, 12326« 

•0,900000 

O' 050000 

9;o12563 

•0,«00000 

0, lOOOOO 

5,999108 

•0,700000 

0,150000 

9,599978 

■0,600000 

0,200000 

3,733028 

•0,900000 

ojsooooo 

2,661237 

•0,200000 

0,900000 

1,998137 

-0,000000 

0,200000 

O'SOOOOO 

1,537096 

0,600000 

ij 191726 

0,900000 

0,700000 

0, 919903 

0,600000 

oJbooooo 

oJqooooo 

0,673855 

0,800000 

0,939792 

0,900000 

0.950000 

0.299315 




CHI 

OELTA.CP 

•0,^90000 

o' 005000 

251993959 

•0.9T5000 

0;012900 

16^190030 

.0,^50000 

0,025000 

11,106818 

•O.POOOOO 

0,050000 

7' 997138 

•0,800000 

0,100000 

9,757178 

.0,700000 

o; J50000 

3,530291 

*0,b00000 

o;2ooooo 

2,793999 

•o,aooooo 

0,300000 

1,930276 

•0,200000 

o;«ooooo 

1,926168 

•0.000000 

0,200000 

0,500000 

1,090151 

0,600000 

0,858607 

o,uooooo 

o;?ooooo 

0,6579ai 

0,600000 

0,800000 

0,966127 

0,800000 

0,^00000 

0,900000 

0,303972 

0.950000 

0.2089U1 

TA • O’, 

X 

,923900 

CHI 

oei.ta*cp 

•0, 900000 

0' 005000 

16' 622696 

•0,075000 

0,012500 

10,200699 

•0,090000 

0,025000 

6,869568 

•0,000000 

0,050000 

9 U09828 

•0,600000 

ojiooooo 

2,611207 

•0,700000 

0,150000 

1,825129 

i;3851l6 

0,929651 

■0,600000 

0,200000 

•0,^00000 

0,500000 

•0,200000 

0,900000 

01689179 

•0,000000 

0-200000 

0,500000 

0 533369 

0,600000 

0,908719 

O.flOOOOO 

0,700000 

0;302620 

0,600000 

0,800000 

0,215505 

0,199678 

0,800000 

0,900000 

0,900000 

0.950000 

0.1057«0 



LOADS OK THIN, LIFTING HINfi 
s8««sssx3aiae3sss8&:zsss:s8 


RECTAKeULAB WING AR o 2 


ID! 

B 

n 

T02 

z 

? 

1D5 

B 

9 

104 

C 

18 

OELTAO 

■ 

4,0000 

EPS 

C 

0,5000 

MACH 

e 

0,0000 

JJ 

s 

191 

JJMAK 

m 

191 

KMAX 2 

u' 

1 

KK? 

m 

5 

PPNfW 

m 

T 

MMNtW 

P 

5 

NROWS A 

u 

0 

CwTVPE 

3 

0 

SwTYPE 

a 

f) 

0C9 

a 

TTFEFFFPPF 

flCAS 

B 

FTFFFFFFFf 

MACH 

fB/aj/BREF 

= o'.oooo 
* 1.0000 

C 6 AR/ 8 PEF 

* l.OOOO 

ASPECT 

Ratio « a.oooo 


1 !• 1 i«73 


rUJTPPT.rOR.N4 



weights 


combination t 

r, 000000 0’ 000000 o', 000000 
combination y 

CL a 2,069057 

INDUCED DRAG a 0,97069<l 

VORTF.X DRAG factor a I 000976 

rolling moment M/ (Q*S*2*RmEF 5 s 0,000000 

right ROOT bending HOmEnt MBR/ (Q*8*8REF ) « 0,527946 

LEFT ROOT BENDING MOMENT MBL/ f 0*3*BfiEF ) a 0,927946 

lift on ETA.OT.O a 1,234519 

^ lift on ETA,LT,0 « 1,234519 

center of pressure of right half. Y/BREF a O' 42T653 

centfR of pressure of Left half, t/brff a •olaBTsii 

pitching moment ABOUT X a 0 AND CENTER OF PRESSURE 
J CM/(QaS*CBaR) X(C,P,)/CBaR 

H .0,512218 0,207465 

23 .0,512238 0,207465 

«T -0,512258 o;2o7465 


PTa 

CL*C/2R 

CL 

CL*C/CCL.*C AVG 

2CMf l/4)C/& 

C,P,fl/4) 

combination 

0,980785 

1 

0^158409 

0,613634 

0' 256632 

o', 060014 

•0,094714 

0,923800 

0, 509436 

1,237745 

0,501306 

0,105865 

•0,083914 

0,851470 

0,446211 

t 704924 

0,722922 

0,123848 

•0,069386 

0,707107 

0j563064 

0,655877 

2,252254 

0,912198 

0,123398 

•0,054789 

0 555570 

2,623506 

11062561 

0,112760 

.0,042961 

0,382685 

0,722500 

2.8R0001 

1,170495 

0,101712 

•0,055194 

0, 195090 

0,762329 

1,049315 

1,215021 

0,094968 

•0,031150 

0,000000 

0.775547 

3,102187 

1.256434 

0.092922 

-0,029954 



COMBlNATtONi 


1 


ETA p 0. 

X 

000000 

CHI 

OILTa-CP 

31*478989 

•0, 990000 

o' 005000 

•0,975000 

0,012500 

19,750305 

13,779659 

•0,950000 

0,025000 

•0,900000 

0,050000 

9,481920 

•0,000000 

0,100000 

6.347512 
4* 899052 

•0,700000 

ojisoooo 

«0|800000 

0^200000 

4,004580 

•0,000000 

0' 300000 

2,897952 

•0,200000 

0,400000 

2,205022 

•0,000000 

o;sooooo 

1,712043 

0,200000 

0,400000 

1,332597 

0,400000 

0,700000 

1,021410 

0,400000 

0,800000 

0,748692 

0 eooooo 

0,900000 

0,480876 

0,900000 

o;9soooo 

0.325402 

ETA » 

X 

382700 

CHI 

DEUt A«CP 

.0,990000 

o' 005000 

29' 991379 

•0,975000 

0,012500 

18,800644 

•0,950000 

0,025000 

13,090129 

•0,900000 

0,050000 

8, 988827 

•0,600000 

0 100000 

5,901541 

•0,700000 

0,150000 

4,590199 

•0,400000 

0^200000 

3,730894 

•0,400000 

01300000 

2,670218 

•0,200000 

o;4ooooo 

2,010678 

•0,000000 

0’ 500000 

1,544572 

0,200000 

0,400000 

1,193743 

0,400000 

0,700000 

0,909155 

0,400000 

0^800000 

0,663239 

0,800000 

0^900000 

0,425072 

0,900000 

0,950000 

0,287423 



FT * * o ; 

,707100 

X 

CHI 

• 0,990000 

oloosooo 

• 0,975000 

o '012900 

• 0,950000 

0,025000 

• 0,900000 

0,050000 

• 0,800000 

0,100000 

■ 0,700000 

0,150000 

• 0,600000 

01200000 

■ 0,400000 

0,300000 

• 0,200000 

0,400000 

• 0, 000000 

01500000 

0,200000 

0,600000 

0,400000 

0 700000 

0,600000 

0,800000 

0,800000 

0,900000 

o’.oooooe 

0 ', 95 oOOO 


FT * B 0 

,923900 

X 

CHI 

• 0,990000 

0 ' 005000 

• 0,975000 

0,012500 

• 0,950000 

0,025000 

• 0, 900000 

01050000 

• 0,600000 

0,100000 

■ 0,700000 

0 150000 

• 0,600000 

0,200000 

■ 0,400000 

0^300000 

■ 0,200000 

01400000 

• 0,000000 

0’ 500000 

01200000 

0^600000 

0,400000 

0,700000 

0,600000 

0,800000 

0,800000 

0,900000 

0,900000 

0 ' 950000 


DELTa^CP 

26;253292 
i ^; ja 7 o«i 
11,276095 
7,586061 
<*,852101 
J , 581182 
2,802904 
1,872190 
1,534875 
0,994303 
01766515 
0,604404 
0,474023 
0.337608 
0'.243238 


ocuta-cp 

17,194107 
10,620198 
71216769 
4 706879 
21817730 
11955975 
1,403632 
0,798974 
01495142 
0,540709 
0,275350 
0,256767 
0,252835 
0,225860 
0.180197 



LOADS ON thin, lifting wing 
•sssasscBssscaxSESssssssess 


rectangular wing AR a 2 


IDI 

fi 

0 

IDS 

» 

? 

103 

e 

R 

loa 

9 

19 

OELTAO 

a 

4,0000 

FP8 


0,5000 

HACH 

m 

0,0000 

JJ 

s 

191 

jjhax 

« 

191 

NMAX2 

9 

5 

kk? 

9 

H 

PPnCw 

9 

5 

MMNEW 

• 

n 

NROwBA 

9 

18 

cwtype 

m 

n 

8WTYPE 

m 

0 

BCS 

fi 

TTFf FfFFFf 

BCAS 

■ 

ftffffffff 

MACH 

(B/a)/BREP 

a O'OOOQ 

a 1.0000 

C8AR/BREF 

a 1,0000 

aspect 

RATIO a a'.OOOO 


ii«n«73 


OUTPUT. roR.NS 



‘^EIGHTS 


COMSINaTION 1 

1,000000 0*. OOOOOO 0' 000000 


COMBINATrOf^ 1 

• 2.060oas 

INDUCED drag ■ 0.970690 

VORTfX DRAG factor a 1,000969 

rolling moment m/ f ) 

9IGhT root bending moment m8R/(Q*S*0R£F) 

MOMENT MRL/rO#S*BprF) 

* 1,239529 

a 1,239524 

OP RIGHT HALF» V/9REF 
OF LpFT half, V/8RFF 

0 and CENTER OF 


left Root bending 
lift on fta.gt.o 
LIFT ON ETA,LT!o 
center of pressure 
center of pressure 
pitching moment about X a 

J CM/roASACBARi 

il oO. 510560 

25 *0*510560 

*0.510568 


u 

a 

a 


a 

a 


pressure 
X(C,P.)/C8AR 
0,206788 
0,206780 
0,206768 


0,000000 

0,527907 

0,527907 


0’,927620 

-0,927620 


eta 

CL*C/2b 

CL 

CL*C/CCL*CAVG 

l/4)C/9 


combination 

1 





01960765 

0,925000 

0,031470 

0,707107 

0.555570 

0,302685 

0,195090 

0,000000 

O; 158293 
0,309171 
0,496098 
0,563105 
01655981 
0,722562 
0,762572 
01775507 

0,635171 

1,236685 

1,769190 

2,252919 

2,625929 

2,890298 

3,099988 

3.102397 

01256445 
0,500875 
0,722623 
0,912262 
1,062727 
11170591 
1,235005 
1 .256995 

0,001057 

0,117555 

0,123910 

0,121039 

0,112390 

0,102298 

0,095699 

0,093229 

•0,1 200 1 8 
•0,095055 
-0,069166 
•0,053735 
•0,092636 
•0,035394 
•0,031382 
•0,030049 



COHSINaTTON t 


ETA a O'OOOOOO 


)( 

CHI 

delta.cp 

•0,990000 

0,005000 

31^495756 

•0,975000 

Ojot2500 

19,760391 

•0,950000 

ojoasooo 

13,766187 

•0,900000 

0,050000 

9 487723 

•0,600000 

0,100000 

6J349156 

•0,700000 

0,150000 

4,699653 

•0,000000 

0' 200000 

4,004546 

• 0 ,()00000 

0,300000 

2,697223 

•0,200000 

0,400000 

2,204043 

•0,000000 

0,500000 

1,711021 

0^200000 

0,600000 

i;33l676 

0,400000 

0,700000 

1,020676 

0,600000 

0,600000 

0,800000 

0,748191 

0,900000 

0,460619 

0,900000 

0,950000 

0.325265 

FTA ■ 0, 

X 

382700 

CMI 

DEUTa.CP 

•0,990000 

0’ 005000 

30' 062637 

•0,975000 

0,012500 

16:842285 

•0,950000 

0,025000 

13,125600 

•0,900000 

0,050000 

9'00l7l5 

•0,600000 

0,100000 

5;984263 

•0,700000 

0,150000 

4,508173 

•0,600000 

0 200000 

5;T26260 

•0,400000 

0,300000 

2,66353? 

•0,200000 

o>ooooo 

2,004232 

•0,000000 

0' 500000 

1,541389 

0,200000 

0 600000 

1,190875 

0,400000 

0,700000 

0,906549 

0,600000 

0,600000 

0,664666 

0,800000 

0,900000 

0,427675 

0,900000 

O' 950000 

0.290082 



FT* ■ 0',70T100 



CHI 

DELTa-CP 

•0,990000 

0^005000 

25,829910 

•0,975000 

0,012500 

16;i13669 

•0,9S000Q 

oj025000 

11,135612 

•0,900000 

0,050000 

71519027 

•0,600000 

0,100000 

6 866690 

•0,700000 

0,150000 

3,601092 

•0,b00000 

o:2ooooo 

2,837557 

•0,000000 

0, JOOOOO 

1,917151 

•0,200000 

0,^00000 

1,376267 

•0 , 000000 

ojsooooo 

1,026673 

0',200000 

0,600000 

01762193 

0,000000 

0,700000 

0,606663 

o;,sooooo 

0,600000 

0,900000 

0,660018 

0 ,600000 

0,317323 

0,900000 

0,950000 

0.226656 

tab o’. 

923900 

CHI 

oelta*cp 

•0,990000 

0',005000 

18*919891 

•0,975000 

0,012500 

11,626663 

*0,950000 

0C025000 

7,627620 

•0,900000 

0,050000 

5,006716 

•0,600000 

o; 100000 

2,867151 

•0,700000 

0,150000 

1,868619 

•0,600000 

0,200000 

11272831 

•0,000000 

OJ300000 

0,619891 

•0,200000 

0,600000 

0 323218 

•0,000000 

0,500000 

0,210693 

0,200000 

0’ 600000 

0,201001 

0,600000 

0,700000 

0 262650 

0,600000 

0 800000 

0,291610 

0,600000 

0,900000 

0,2969?7 

0>00000 

0.950000 

01265205 


rograw Halted due to stop command 
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onnoon nnr^ 


integer UC I ,UC0»U6»U7,Ul2,CWTYPEr SWTYPE,UNSYM, T JJH aX, 

I PPNfcW,PP,OoISK 

real MaCH 

OOUhLE precision OUMI, DUM2, ET*, STHETa 
LOGICAL BATCH,CONV,BC8,BCaS,8YM,REa07,OPEN<,, CHECK 

b I HE’NS 1 ON T I TIE ( 26 ) # BNK ( 630 ) # »NK 1 { 750 ) » CH ICP f 20 ) » 

IBCSf l0)»BCASM0)»NlNOEX(<l7)#S7ORl(3fl3),SToRpc383)# 

2STHETAf3a«),X8ILIP(38l)#COROIPf383)»ETAOc^0)»WpIGHT(50)# 

3ETA(383> 

...NBNKl and NHEI are the dimensions OF {BN«l) AND (HEIGHT), 

DATA UCI ,UC0,U7,Ui2rNBNK 1 ,NWEI ,COM/5»b#7, I2,750f 5o, 3H / 

establish whether THIS IS A BATCH OR CONVERSATIONAL JOB 
WRITE CUCO»8o2) 

REA0(UC1#803) batch 
COnV c , NOT, batch 
lOlR * -PRO 

OPENfe INDICATES whether OR NOT UNIT U6 IS OPfN‘, 

THIS Variable is not needed in standard fortran, 

OPEN6 ■ .EAlSE, 
to CONTINUE 

call 0BEV(t6^t6HRELEASE PTo«EOOj ) 

IF(CONV) WRlTg (UCO,805) 

ODISK determines MHETHfR OR NOT TO PUT OUTPUT ON oISK FILE (UNIT 4) 
IF OpISKeO, OUTPUT WILL pE ON UNIT 6 

IF OdISK ,nE, 0 OuTPijT WILL rE ON tjNlT « ON a FILE wHOSE NamE 
IS generated by THIS program; 

REa 0(UCW810) ODISK 
OOI5K«MOO(ODI5K, to) 

ENTER ODISK, LT,0 TO TERMINATE EXECUTION, 

if(Odisk.lt,0) stop 

U6 * UCO 



nr»r%r»nr* rarionrin n n ooorto 


IP(0DISk ,NE,0)U6«« 

0PIN6 ■ UC0.NE.U6 

IF f .NOT.OPpNf,) GO TO 2o 

THE CODE below JS FOR THE a«ES ' TSS VERSION^ OrEY GIVES COMMANDS 

, ,,,,T0 THE TSS operating SYSTEM, THE VaRIaBLE OOISK IS NOT 

,,,, .necessary on a standard version, 

WRITE(UC0#8U) ODISK 
CALL CVHT COOISKi 1 , 

I «aH(»D0FF FT0aF00I», OUTPUT, pOR.Ni ,11, AJO 
ISTORI >8^SH(6Att} ) 

Call obey r32»ST0Ri) 


20 JFfCONV) WRlTEfUcO,8l6) 

RCADcUCIfSlo) IDDU10D2#I003#10D<1 

check « 1001, GT.O 

read? * iOiR.NE, iDOl 

IDl « lODl 

ID? 9 ID02 

10) a IDD) 

IO(i a IDDti 

IF ( ,NOT ,«EaD7) go to «fl 

OPEN and read first RECORD OF THE GEOMETRY FILE, 

...FOR aMES» TSS VERSION ON|,Y, GFMFIL ISSUES OOEF COMMANDS TO THE 
...TSS operating SYSTEM, 

Call GEMFlLdOn 

REAO(UT) 10#PP,MM,CWTYP|,8WTYPE,UNSVM,NOL,NOT,MREP, jjmax, 
INFLaPS.T nifc.NTITL 

NDu X number of DISCONTINUITIES ON LEADING EDGE 
(MINIMUM Value of n 

NDT ■ NUM8FR OF DISCONTINUITIES ON TRAILING EDGE 

(MINIMUM Value of n 

PP a original number OF CHORDWISE CONTROL POINTS 
MM 8 original number of SPaNWISE CONTROL POINTS ON 



rifioonn r»nnr* 


f entire wing. 

c 

C CHECK ID number 

IF ( .NOT, CHECK) GO TO SO 
IF(ID1 ,EQ,IO) GO TO 30 
1 E(BaTCMJ C*Ll 8 T 0 P 2 (UC 0 # 

p 1001 DOES NOT Match geom pile Value *, 

jFLOATdOn 

PAUSE «1001 DOES NOT MaTcH GEOM FILE VaLUE* 

30 CONTINUE 
lOtR • IDI 
C 

c OPEN and read first record of solution file 

«0 CONTINUE 

*FOR aMESi TSS version ONLY, rNkFIL GIVES OOEF COMMaNOS TO THE TSS 
.OPERATING SYSTEM, 

CALL BNKFILfl 0 i»I 02 dO 3 »lO<i) 

REA 0 {Ut 2 ) lOl , 102 , 103 # 104 , title, SYM. NS YM.NaSYM, 

J(0C8(I),I»1, lO), feCAS(I),l«j,iO),NMAX2,KKg,JJ, 
2PPNEH,mMNEW,NR0WSa,0ElTA,EP5,hACH 

NMAX? « NEW number OF C.H 0 R 0 WISE MODES 
KK^ • NEW number of SPaNWISE MODES 
PPNEW • NEW number of CHOROWISE CONTROL POINTS 
MMNEW p NEW number OF SPANWISE CONTROL POINTS 

IF f. NOT, CHECK) GO TO 60 
C 

c identification checking 
IF flOOl . EG, Ion GO TO 50 
IF(bATCH) call STOPEfUcO, 
n looi DOES NOT Match solution file value*, 

SFLOATdOd) 

PAUSE 'IDDI does NOT MaTj-H SOLUTION FILE VALUE’ 

50 IFd 002 ,e 0 .l 02 } GO To 52 
IFfBATCH) CALL ST 0 P 2 fUC 0 , 

1' ID02 oOfS NOT Match solution file Value*, 
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) 


2FL04T(1D2)) 

Pause fioDj does not match solution file value* 

5? IF(1DD5,E0,ID3) go TO 5i| 

IF(BaTCh) call STOPSfUCO, 
p IDDJ DOES NOT MaTCH SOLUTION FILE VaLUE», 

2FL0aT{ 10}} ) 

Pause 'Iddt does not maTch solution file value> 

50 1F( lODO.eo, l0«) Go Tn 60 

IFfBATCH) CALL STOP^tUCO, 

P 1000 DOES NOT Hatch solution file Value*, 

2FL0aT(I00} ) 

Pause *ioo« does not maTch solution pile value* 

60 CONTINUE 

r. 

C PRINT HEaOING 
WRITE (U6rB"t9) 

WRITF(U6,620)(TITLE<n»I«l,NTlTL)f 101 #102,105^0^# 
1 DELTa,EPS,HaCH,JJ,JJMaX,NMaX2,(<K2,PPNEW,MHNEW,NR0WSa,CWTYPE, 
^ ?SWTYPE, (6CS(n#l*l»lO), (BCAS(n,I#l,10) 

(O £ 

NCDIM , NSYM ♦ NASYH 
NDLPT ■ NOL * NOT 

,LLHAXI and LLMAX? A«E two LIHITS for the number of 
.combinations in subroutine SBLOaD, 

LLMAX2 c NBNkI/ (NMAm2 * KK2) 

UMAXI s NwEI/ NcOIM 
TJJHAX * 2 * UJHAX 

CALL SBlOAD(UcPU6,UcO,U7,U12,5YM,NSYM,NASYM,NpoIM, 
1 NMaX2,KK2, JJMaX,PP,MM,MREF ,NDLPT,TJJMaX,«EAD7, 

20NK,bNKi ,CHIcP,NINO£x,STO«|,ST0R2,ETA,3TH£Ta, 
jXSlLrP,COROIP,FTAO,CONv,WFlGMT,LLMAXl,LLMAX2,COM,MACHl 
GO TO 10 
C 

802 FORMaTUH PENTER BATCH* } 

80 l FORMATfUl) 

Bo5 fORmatuh penter ooism (neg. halts )' ) 

810 F0RMATfl6I5} 

Rl2 FORMATrgOH OUTPUT IS ON 


OUTPUT, for, N pn 
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(i\h FORMaTHH , 'ENTER lOlf 102, IOS,lDt)* ) 

619 FORMATC MLOADS ONi THIN, LIFTING WIng?/ 

1 > •aBamt8CBiiBsaBiBB«iss>t«aca • /) 

820 FORMATf lx,20A«//,i 101 a ID2 s f, 

U}/,« 103 a ',13/, • 104 s *,I3/, 

DeLTAO b •,F8,4/#' fFS a ',F6,4/,( HACH a 

3F8,4/,r JJ a ',13/' JJMAX b »,I5/* NMAX2 s ',I3/, 

4» KK? e PPnEW a ',15/, • MMNEW c 

5l3/,’ NRORSA » »,i3/,* CWTYPE a 

6' SWTVPE a ',13/, ♦ BCS a »,10Ll/» BCAS a ',10LD 

C 

End 

subroutine S8t0AD(UCI,U6,UC0,U7,Ul2,SYM,N8YM,NASYM,NC0IH, 
1 NMaY,KK, JJMaY,PP,MM,MREF,N0LPT,TJJMaX,«£A07,BNK, 
2rNk1,CHIcP,NIN 0EX,5T0R1 ,ST0R2,eTA,8TM£TA,X8lLlP,C0RDIP, 
3 £TaD, CONV, WEIGHT, LIH 4X1 ,LIMaX 2 , COM, MaCH) 

C 

double precision ETA, STHETa,PID,TNETSW,OTHETA, theta, 
ITH£TcW,THET 8 c,DELTHE 
integer U6,UrO,U7,UU,PP,TJJMAX,UCl 

lOClCAL 5YM,CONV,PRtSUR,PHINT,SLOAO, ASL0 a0,TE8T,T£STN1, 

1 TbSTN?, RE A07,C0NT IN, half, TEST?, TESTS 
REAL MACH 

real NET$, news, KNEs 

c 

DIMENSION BNK(NMaX,KK,NCOIM),BNKI (NM aX , KK , LLHAXg ) , 
1WEIGHT(NC0 IM,LLMaX, ),CH lCP(PP),NlND£XfMM),STORj (JJM aX), 
2ST0R2(JJMAXJ ,t TA(,JJMaX), 8 ThpTa< JJMAX),X 8 ILIP(JJMAX), 
SETaO(NOLPT) ,THETSW f200 ) ,THEtcW( lOO),THETSC(475,CAROIMe«) , 
aSLOAOno) ,asloaO{io) ,coRoiP( JJMaX) 

N0IMl,N0IM2 aNO NOIMJ aRE THE DIMENSIONS OF THfTSW, 
thetcw, and THETSC 


DIMENSION TSScOM(2o) 

data N01M!,NOIM2,N01M3/200,100,47/ 

DATA PlD,Pl2,F0Pl,F0Pl,PlS.Pia/S.l«l592S53589793 DO, 
11,570796, 1,271240,2,546479,0,5926991,0,7853982/ 

Data Etas / iheta / 



data XCPt / JHXCp / 
data spa* / 5HSPA / 

DAT* tCP* f 3HECP / 

data net* / SHNET / 

data PRi* / 3HpRC / 

OATA_Fl.A* ( 5HFj,A / 

Data con* / 3hcon / 
data STOS / 3HST0 / 

DATA new* / 3HnFw / 

DATA MEl* / 3HWEI / 

DATA KNE* / *HKN£ / 
data PWE* / 3hPWF / 
data TSS* / 3HTSS / 

C 

C read second record of geometry FIUE 

ir>REA07)REA0(U7) CHI cP » Nl NOEX , STOR I , STOR?, STQR 1 , SToRj , 
JETA,STHETA, 

IXSILIP,CONDIP,0RaT1O,C8aRBR, AR,TR,ZmaCH,ETaD 

o» C 

^ c WRITE balance of heading 

write rUb# 6000} MACH^BRATlOtCBARBR , an 
C 
C 

CONTIN a cony 
3 CONTINUE 

IF(SVM) go to 3 
C 

c reading solutions for UNSYMMETRIc wing 

DO gOOl I»l,NCOIM 

2001 REaD{U12) ( (8Nk<N,k# n »NaUNMAX) ,k«WKk) 

GO TO 11 

5 IF(N8YM,fO,0) go to b 
C 

c reading symmetric solutions for Symmetric wing 

DO gOOg IslfNSYM 

2002 read (U12) C(BNK{N»K» l},N3l,NMAX),K»lfKKr2) 

DO a 1 8 l.NSYM 

DO 4 K a 2fKK|g 
pO 4 N a 1 ,NmAx 
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c 

c 


<1 BNK(N,K, I) « 0« 

6 IF fMASYM.EO.O J GO TO Jt 

U • NSVM ♦ 1 

reading antisymmetric solutions for symmetric NIMG 

00 a005 Il21,NCD|M 
2003 REaO(UU) ((BNK(N,K/n»N»l,NMAx)fK«2»KK#2l 
00 7 1 s n,NCDIM 
DO 7 K « 1,KK,2 
00 7 N • l,NwAX 

7 8NK(N,K,n « 0, 

\\ continue 

C 

TfSTNl s NmAX,GT, 1 
TESTNg c NMAX,GT*2 
■jUf^PR a fJjMAX^n / 2 
OELTH£aPIO/OFLOAT(JjMAX*l ) 

ARFFF a AR / BRATI0**2 
ot C6ARe^ • CSARBR * BRaTiO 

*" CONI a PI2 * ARIFF 

C0N2 a CONI / 2, 

CON3 • C0N2/2^ 

JJHAxi » JJMax * 1 
ISUB a JJMAX) * JJMAX} 

the dimension of stmeta in the Calling program must 
BE ONE Larger than jjmax 
STHETA{JJMAX| ) • 0,00 

enter combination code 
defines a new set of weights 

15 IF(CONV) WRITE(UC0,BM5) 

LCOmb.GT.O, then LCOMBiTHE NUMBER OF COMBINATIONS AND 

weights of Each combination will be entered by the 

. , , , .USER, 

IF LCcoM,Eo,o, Then the number of combinations a The number Of 

Cases (Solutionsj and each solution with a factor of i o 
,,,,,15 treated as a combination. 
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IF LCOMB,tO,-l, THEN THE NUMBE« OP COMBINATIONS « THE NUMBER 

symmetric Cases anO each symmetric solution with a 
,,,,,FaCTOR of 1,0 IS treated as a combination, 

LC0Mb.E0,,2. then the above holds, except THaT the 

antisymmetric CaSeS are involved; 

UCOMg,LE,-S, THEN the CURRENT SET OP WEIGHTS WILL BE 

♦••tfUSED, 

REAOfUCl*e020) LCOMH 
IP(lABSaC0HB4>n ,LT,2) go to 32 
IF(LcOM8.GT,0)U.«LrOMB 

ife if(LL.le,mino(LLmaxi ,llmax2)) go to it 
IF (,NOT,CONV) Call 5T0P?(UC0, 

11 TOO Many combinations »»FLOATai)) 

Pause * too many combinations • 
r*o to is 

17 IFrLCOMB.Lt.-3) GO TO 3’ 

It s NSYMaI 

o> DO jO L » IfLL 

IFCCONV) WHITCrUcO.8025) t 
IF fNSVM.rO.O) GO TO 20 

iFrcONv) wRiTErUro.eoiO) nsvm 

RE AOfUCl rS035J (WEIGHT (I#L)« I*t#NSVM) 

20 CONTINUE 

1E(NaSYM,EO,{)) GO TO 30 
IF(CONV) WHITEfUcO»eO<<0) nasym 

REAOf UC 1.80 3B) (WEIGHT r I,L) . Isll .NCOIM) 

30 continue 
GO TO 59 
C 

32 IFacOMB.EO.O) U B NCOIM 
rF(LCOMB,EQ,-i) LL • NSYM 
IF(LC0 Mr,EQ,.2) LL « NaSYM 
If(CONv) GO TO 33 

if(LL,gt,hino(LLMaxi,llmaX2)) Call stop?(uco. 

V TOO many cases* .FLOaT(LL) ) 

IF(LL,LE,0) CALL STOP2(UC0, 

1' number OF Cases is not positive i.floaT(lU) 

GO TO 5^ 
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33 rONTlNUE 

lF(LL»GT,MlNO(tl„MAXl ,tLMAX2) 1 
IPAUSE tTOO MANY CASES» 

IPeLL.Lfc.O) Pause 'number of cases is not POSITIVfi 
5« continue 

35 t «_ l,u 

DO 35 I s tfNCDIM 
55 WCIGHT(I,L) 8 0« 

1 » 0 

iFat.eo,-?) I It NSYM 

DO 36 L«1,LL 

r e r + 1 

56 WElGHTfl,!.) c i„ 

59 CONTINUE 

IF fCONV,ANO,UCO,,£Q,U6) 00 TO tt9 
SET OF WEIGHTS IS TO «E PRINTED EXrEPT IN 

».,,,the conversatiomau mode ano when all 

♦.♦..OUTPUT IS ON ThE terminal 

wRITF(U6,80«5) 

DO 65 L«1.LL 

65 NRITfe (U6.90S1)L. ( Mr IGMT ( I . L ) ♦ I« I . NCO lM> 

69 CONTINUE 

NOW tHaT the WEIGHTS HaVE BEEN DETERMINED. THE 

.....COEFFICIENTS OF EaCH COMBINATION WILL BE CALCULATED 
DO 150 L s l.LL * 

00 160 K B l.KK 
00 160 N 8 I.NmAx 
160 BNKifN,«,L) 8 0, 

00 150 I 8 l.NCOfM 
OUMl • WElGHTfl.U 
00 150 K 8 I.KK 
DO ISO N 8 I.NhAk 

150 BNK1(N,K,L) a 0NKlfN,K,L) ♦ OUM} * BNKfN.K.I) 

TEST 8 KK,LT,2 
00 190 L ■ l.LL 
SLOADa) ■ TrST 
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JF(SuOaD(L)) go to 175 
SLOAOa) s BNKl(l,2,L),eQ.O, 

1F( ,NOT,SLOaO(U 1 GO TO i 75 
DO 170 K • 2|KK,2 

DO 170 N » I.NhAx 

170 SlOAOtU • 8lOAOfUJ,ANO,8NKl(N,K#t),Ea,0. 

175 ASLOADfU « .NOT , SLOAO f L ) 

IF(8L0aD(U) go to 165 
A$LOAD(U • eNKl{l,l,U,eQ.O, 

IFt,N0T,A5l0A0(l>) GO TO 185 
00 160 K ■ 1 |KK ,2 

00 160 N s I.NmAx 

160 ASUOAD(L) e ASLOAO f 1) . AND, BNK 1 f N, K , L ) , EO* 0 , 

185 continue 

SIOaO(U nILL be TwUE if THE WING IS SYMMETRICALLY LOADED, 
asloaO(L) will be true if the wing 18 antisymmetrically 
loaded. 

190 CONTINUE 

IF (rOM ^EO^ WEI*) 60 TO 200 
IF (rOM ,EQ, KN£S) GO TO 200 


f THET 3 W) , f THETCW) 4 NO (ThETCS) a«E THE aNGULaR POSITIONS 
AT which SPaNWISE quantities and the CHORDWISE and 
SPaNWISE POSITIONS at WHICH PRESSURES WRL rE COMPUTED 
respectively; these arrays can be changed by THE 
»etaS», 'xcp’# and >ecp* commands respectively. These 
ARRAYS are ONLY USED IF a tSPANLOAOt OR 'PRESSURE* 

command IS Later given. 

These arrays are derived from thE original spanwise 

CONTROL POINTS ON THE GEOMTERY FILE, 

IF(MM,GT,NDIM 1 ) CALL STOPgfUb. 
l» NOT FN0U6H ROOM FOR DEFAULT SPANWISE LOAD STATION. 
^increase NDlMl TO 'iFLOAT(MM)) 

IF eMM,GT,N 0 lM 2 ) CALL STOP2fUCO, 

1*N0T ROOM FOR default CHORDWISE PRESSURE STATIONS, 

2INi“, NDIH2 TO »,FLOaT(MM)) 



ir(MM,GT,NDlM5) CaU STOP^fUCO# 

I 'NOT ROOM FOR 0£FaUlT SPANWISE PRESSURE STATIONS, 
2INC, NOJMJ TO '#F£OaT(MM)) 

JRATIO a JJMAXl / fMREFtJ) 

OO 195 H ■ 1,MM 

index ■ N1NDEJ({M:i * JRATIO 

THETA a OFIOAT ( INDEX ) ADEL THE 
TKfTSW(M) a THfTA 
ThETCW(M) a THETA 
THETSCfH) a THETA 
195 CONTINUE 

NTASS a HH 
NTA6C a MH 
NTA6SC a HM 
50 continue 

200 continue 

IFrcONV)WRITE(UCOiBOlO) 

C ENTER COMMANDS 

REAOfUCIf eo<f> 0 ) (CAROlMd) , lai ,4) 

COM a CARDjM(l) 


IF 

(COM 

.£Q. 

ETah) 

CO 

TO 

500 

IF 

(COM 

.to. 

XCP$) 

GO 

TO 

522 

IF 

(COM 

. EO » 

SPaSI) 

GO 

TO 

542 

IF 

(COM 

« E Q . 

EcPil) 

Co 

To 

56ft 

IF 

(COM 

.to. 

NETji) 

CO 

TO 

380 

IF 

(COM 

.to. 

PHEJi) 

GO 

TO 

400 

IF 

(COM 

.EO. 

FLAIi) 

GO 

TO 

440 

IF 

(COM 

.EO. 

CONf.) 

GO 

TO 

460 

IF 

(COM 

.EG, 

STO$) 

GO 

TO 

46o 

IF 

(COM 

.EG. 

NEWf) 

GO 

TO 

485 

IF 

(COM 

.Eg. 

WElf,} 

GO 

TO 

15 

IF 

(COM 

.tc. 

KNE5.) 

GO 

TO 

490 

IF 

(COM 

.EG. 

PWE*) 

GO 

TO 

500 

IF 

(COM 

«EQ. 

TS5!,) 

00 

TO 

520 


WRITE (UCO,B065) ( C ARDIM ( I ) # I * 1 » 4 1 

IF(cOntIN) 00 TO 200 

STOP 


C 
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rr a8-*-*6IVES ETaS for printing the LOaOS, 

THE RESUUTS WRU BE TME Sa^E AS IF 
NTyPE s 0 

ENTER Value for ntype 

100 IF(CONV) WRITE(UCOf 8070) 

RtADfUCl*8020) NTYPE 
IF fNTYPE.EO.O) GO TO 30l 
IFfNTYPf.EO.n GO TO 50 « 

IFfNTYPE.EO.?) 60 TO 508 
IF{NTYPE,E0,5) go to 5l^ 

IF (NTYPE, EG, 5) GO TO ll7 
IF(,NOT,CONV) Call 5T0P2(UC0, 
p NTYPE invalid «#FLOaT(NTYPE)) 
NRITE(UCO#80TBi) 

GO TO 300 
C 

NTYPE I rt 

^ 301 NO • ( JJMAX+1 ) /(MReF^I ) 

GO To 3fl5 
C 

C NTYPE 8 1 

SOU IF(CONV) WRITE(Ur.O,8080) 

WEA0(UCI|80^0) NO 
505 n *0 

U 9 (JJMAy ♦ 1) / no 
NTABS c T2 • i 
C IF(HALf) NTABS » 12/2 

DO 307 J « NO,JUPPR,NO 

11 ■ n ♦ I 

THfTSN(Il) ■ dFLOAT(J)*dELTHE 
C IF(HaLF) go to 307 

12 e 12 • 1 

THETSW(I2) a P I D-1 HE TSw ( 1 1) 

3o7 CONTINUE 
GO TO 200 
C 

C NTYPE 8 ? 


IF NOT GIVEN 

given with 



30P IPfCONV) WHITErUCOf 8O0SJ 
«EAD(:UCI,fio20) NSTa 
12 s NSlfA ♦ I 
NUPPR ■ i2/2 
NTABS c NSTA 

C IP (HALF]) NTABS m NUPPR 

DtHETA ,1 PlD / DFLOATfNSTA+n 
DO 310 n s l,NUPPR 
THeTSW(I|) s DFLOAT(H) * OTHt'TA 
C IF(HaUF) go to 3110 

12 • 12 . 1 

THETSWflpj s PIO • ThET5W(Ii) 

31ft continue: 

CO TO 200 

c 

C NTyPE ■ % 

312 nTabS m ft 

C ENTER TAStE OP VALUES FOR ETa 

IF(CONV) NRIT|(UC0,8090) 

313 NTABS B nTABS t I 
READfUCI»803S) TheT5W(NT A0S) 

IF (TH£TSW(NTa 8S) „LT*1 tOO) GO TO 3l3 
IF (NTaBS,LE,NOIM.| ) GO TO 3i4 
CALL ST0P2fUC0» 

1' TOO many entries in spanuoao eta Table 

2 SNGLfTHETSWfNTABIS;n) 

3H nTAbS 8 NTAbS • 1 
DO 316 N 8 I.nTABS 

316 THETSWfN) e OARCOS f T h£ TS« (N J ) 

GO TO 200 

C 

C NTYPE ■ % 

317 IFfCONV) WRITEfUcO,8087) 

HEaOcUCi » 803S) £TM|N,ETMAX,0ETA 
IF(DETA,EO,0) OETA s ETMJN 

IF (ETMaX ,EO, 0) ETMA* « ,9999 
thpTSNn) B pTMlM 
NTABS 8 1 

3ift IFfTHETSN(NTABS) ,GT,ET maX) GO TO 319 
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NI.AST z NTAbS 
NTABS s mTAHS ^ \ 

THgTSW(NTABS) « THfTSWfNLAST) + dETA 
CO TO 5l8 

319 1F(NTABS,GT NOIMDCALL 8T0P?(UC0# 

jiSPANWlSf load TABUe SIZE e^<CEEDEO », FLOAT < NT A8S ) ) 
NTABS ■ nTABS • 1 
DO 321 I « 1. NTABS 
321 THETSWfl) s OARCOS(Th£TSW(I)) 

CO TO 200 


XCP-.OEF INES KtS FOR PRINTING THE fHORDWiSE PRESSURE 
DISTRIBUTIONS 

32? IF(CONV) WRITE(UrO,8oTo) 

«EAO(UCI»«020) NTYPe 
IF fNTYPE.tO.O) GO TO 5?3 
IFtNTYPEjEO.l) GO TO 3?« 

IF (NTYPE.EQ,?) gO TO 3?S 
1F(NTYPF.E0.5) go to 332 
1F(NTYPE,EQ««) go to 332 
lFfNTYPE,E0.5) GO TO 338 
IF t ,NOT.CONVi Cali 8T0P2(UC0» 

1» ntYPE invalid wFlOaTCNTYPET) 
wPITf. fUCO»8075) 

GO TO 322 
C 

C NTYPE ■ rt 

323 NO ■ { JJMAX + n /(MR fF ♦ U 
CO To 3?5 
C 

C NTYPE s 1 

3?a IF(cONV) WRlTEfUcO,e080) 

Rf AD(UC i . 8020 ) NO 
325 n *0 

12 z (JJMAY 1) / no 

NTABC m 12 • 1 

00 327 j * NO,jUPPP,NO 

n R n ♦ 1 



THEtCW^ID tt DrLO*T(J)*OELTHE 

U * U . I 

THETCWfl?) n PID - THETCW(IJ) 

327 CONTINUE 
GO TO 200 
C 

C NTYPE a ? 

32B IFfCONV) WRITEfUc0#808«^) 

RE AD rUCi ^ 6020) N$Ta 
12 a NSTA 4 1 
NUPPR a 12/2 
NTA6C a NSTA 

OTHETA a P10/DFL0AT(NSTA ♦ 1) 

00 5J0 11 a I, NUPPR 
THETCW(ll) a OPIOAT(II) * OTHETa 
12 a 12 . 1 

THETCt^rl?) a PIO » THETCW(It) 

330 CONTINUE 
N4 GOT0200 

® C 

C NTVPE a % OP « 

332 NTABC a o 

IP rC0NV,AN0,NTVPe„E(J.5) WPlTE(UcO,809n 
IP fCONV,ANO,NTVPE„EO.a) WRI TE (UCO, 8092) 
3 J(J NTABCaNTABC + l 

R£AO(UCI,e03S) THILTCW(nTaBC) 
IF(THETCN(NTaBC).LT,J ,) GO TO 33« 
IPfNTABC.GT.NOIMj) CaLL STOP^iUCO, 
i» TOO Many entrjES in t or ChI Table*, 
2SNGL(THETCW(NTa8C))) 

NTABC a NTABC -1 
iFfNTYPE.EQ.U) fiO TO 536 
no 335 N a 1,nTA0C 

335 THETCN^N) a OARC OS (.THE Tc^ f N) ) 
f,0 TO 200 

336 00 337 N • I ,nTA8C 

337 THETCW(N) a OARCOS ( 1 ,00-2, DO * THETCW(N)) 
GO TO 200 

C 



c ntype » s 

35B If(CONV) WRIT6(UC0,80<>5) 

RE*DfOCI,8oi5) XMIN,xMaX,DX 
If fXMAX.EOfO,) XMAX • ,9999 
If(OXtEU,0,) DX«1,*XM1N 
NTA 6 C « 1 

THETCWfNT*BC5 « XMIN 
33R IF(TH£TCW(NTabC) ,(»T|XMaX5 gO TO 3ao 
NIAST • NTABC 
NtABC 8 NTABC ^ 1 

THeTCW(NTABC) n THfTCW f NL AST ) ♦ !)t 
GO TO S59 

3ao If (NT a8C,GT,NDIH2) CaLL STOP2fUCO^ 

PCHOROwISe pressure TaBIe SI2e exceeded • #FLOaT(NTa8C) ) 
NTABC I nTABC • \ 

DO 3ai I t 1, NTABC 
3ai THETCRfn* DARCOS(-THeTCW(I)) 

GO TO 200 

c 

C SPANiOaDS--«PRINT5 results aT SPANWISE stations 
c defined by ETaS COHMaND'. if the ETaS COHMaNO has 

c NOT PEEN given TmE EFFECT IS ThI SAME AS iF iT WERE 

C GIVEN RITH NTYPE s 0 

3a? wRITF (U br^OOO) 

DO 356 I m 1,LL 
WRITE (U6,90Ia) L 
CLTftMsBNKl f 1 , I »L) 

IF(TeSTNi J CLTHM,CLTRM+,5*bNk 1 f2« 

TEST3«CLTHM,NE,0. 

CLCCLC«-l,e-10 

TEST t SlOAOfL) .OR.ASLOADfL) 

KL a i 
KJUHP « 1 

iF(ASLOAoan Ki » ? 

IF(TFST) KJUMP K ? 

DO 355 N « I,NTABS 
THETA a THETSWfN) 

IF fTEST.ANn, THETA, GT, (PI2+I ,n*«) ) GO TO 355 
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c QmOMNT s OUARTtR CHORD MOMENT 

OMOMNT s 0, 

GAMHAtO, 

00 S«5 K a KL,KK#KjUHP 
SKTHET a OSINCDFLOaT(K)*THETA) 
gamma a gamma 4 $KTHpT*6NK 1 { t » K , L) 
lF ('TF8tNl ) GAMMA a GAMMA ♦ SKTHET* , «S*BNK i (2, K, L ) 
IF(TfSTNi)OM 0MNT a QMOMNT S*SKTMET*rNk j ( 2 , k , L ) 

IF (TfSTN 2)0M0MNT a OmOMNT ♦ , 5*SKTHf TtBNK 1 ( 3 , K , t ) 

3«s continue 

GAMMA a NOND I MENSil ONAL CIRCULATION a (CL*C)/f2P) 

OMOMNT a NONDIMEMSIONAL PITCHING MOMENT LOaO 
8{2*CM(1/6i)*C)/b 

ETAS a DCOSfTHiTA) 

A uneaR Interpolation will be used to obtain the 
CH0R0/(B/2) aT the current ETaS station 

(ETa) is in descending order 

IPfETAS,LT,ETAf2)«l ,0*6) GO TO SRR 
s 2 

CO TO 35? 

SaQ CONTINUE 

00 350 J * 3,JJMAx 
J1 a J 

IF(ETAS.GT,ETA(J)) go TO 352 
550 CONTINUE 
55? CONTINUE 

CORDB? a CORDIP(Ji) ♦ f COROIP( J j ) -CORDIP f J 1 - I ) ) / 

I fETA(jn-ETATjl*l))*(ETAS-eT4{ jU) 

C 

C SECTIONAL LIFT COEFFICIENT 

CL a «,*GAMMA/C0R0B2 
C 

c overall lift coefficient 
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aspect Ratio from GEOM file is a(BREF**2)/AHtA, 

CCL a fPI/2)*(B**2/AREA)*(BNKl n,lfU + ,5RNKl(2,l,L)). 

BRaTIO from GeOH file is BREF/fB/^), 

THUS 

(CL*C)/fCCL*CAVG)anAMMA*a/Rl*B«ATlO/tBHKt n, + (2, 1 ) ) 

FOPI a R'./Pl 

IF fTESTS) CLCCLCaGAMMA*F0PI*8RATl0/CLT»M 

CENTER OF PRESSURE RELATIVE TO THE LOCAL QUARTER CHORD 
TEST? a ABSCGAMMA) ,CT, ( ,t-5*ABSfQM0MNT) 
iFfTEST?) COP a-, 25* QMOMNT /gAMMa 

IF(TeST2) WWITF. (U6,9005) E T AS , GAMMa , CL , CLCClC , QMOMNT , 

ICOP 

IF(,N0T,TEST2)WR1TF(U6,9005) ET AS,CAMMa,CL#CLCCLC, 

J QMOHNT 

U55 continue 
S56 CONTINUE 
GO TO 200 


ECP— -define ETaS for printing THE CHOROWISE PRESSURE 

360 IFfCONV) WRITEfUcO.BoTo) 

REAO(UCIi8020) NTYPe 
IF (NTYPE.FO.O) go TO ibl 
IF fNTYPE.EO, 1 ) GO TO 
IF(nTYPE,E0,2) Go To 368 
IF fNTYPE.EO, 3J GO TO 372 
IF(NTYPE,E0,5) CO TO 376 
IF(CONV) NRITEfUC0,807B) 

IFfCONV) GO TO 560 

Call ST0P2TUC0.t NTYPE INVaL ID » , FlOaT (NTYPE ) ) 

C 

C NTVPt s 0 

361 NO a (JJMAx+l) / fMRFF+1) 

GO TO 3«>5 

C 

C NTYPE a 1 

36U IFfcONV) WRITEfUcOiSOSO) 



HEADfUCI f A020) N(3 
365 II • 0 

I? m (JjMAX4t)/N() 

NTABSC R 12 • t 
C lFfHALF)NTAB$C • 12/2 

DO 367 J « NQ,JUPPR|NQ 
It • It ^ 1 

THeTSCflU c DFLOATf J)*neLTHE 
C IF (HALF) GO TO 367 

12 a 12 « 1 

TMETSC(I2J a PID»ThET8c ( I t) 

367 CONTINUE 
GO TO 200 
C 

C NTVPE a ^ 

566 IF fCbNV)WRlTE(UCb#8085) 
REAO(UCJi 6020)NSTA 
12 • NSTA ♦ 1 
NUPPR a 12/2 
NTA6SC a NSTA 

C IFIHAIF) NTAeSC « NUPPH 

dTHETA * PlD / OFLOAT(N8TA4|) 

DO 370 It a l,NUPPP 
THETSCfli) a DFlOATdl) * OTHET* 

C iFfHAlF) CO TO 370 

12 a 12 . 1 

TM6TSC(I2) a PID • TMETSCtll) 

570 CONTINUE 
GO TO 200 
C 

C NTvPE a 3 

372 NTABSC a 0 

c ' ENTER Table of values for eta 

IF(CONV) WRITEcUcO,8090) 

173 NTABSC a NTABSC + \ 

REAoeUCI,e035) THrTSC(NTABSC) 
IF(THETSCtNTABSC),LT,l* J GO To s73 

IF(NTaB 5C,L£.N01M3) GO TO 37a 
IF( ,NOT,cONV) CALL ST0P2fUc0, 
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1 ' TOO many Entries in pressure ETa TABtE'» 
aSNGl{THET5C(NTABSCn ) 

Pause tTou m^ny entries in pressure eta table' 
nu ntabsc s ntabsc » i 
DO 375 N > I, ntabsc 

575 THETSCtN) m OARCOS (THETSC (N) ) 

GO to 200 

C 

C NTYPE « «5 

576 IF(CONV) WRI7EfUcO»80STJ 
READ{UCI»6035) ETHiN,ETMaX,OETA 
IFfO£TA,EO,0,) OETA « ETmIN 

IF f ETMAX ,EQ, 0 , 1 ETmAY s '9999 
TMETSCM) e ETHIN 
NTABSC « I 

377 IF (THETSC(NTaBSC) ,GT,ETMaX) GO TO 578 
NL*5T n ntabsc 
NTABSC 8 NTABSC 4 1 
THETsCfNTABSC) i ThETSC ( NL AST ) 4 OETA 
GO TO 577 

578 IF (NTABSC, GT,N0IM3^ CALL STOp2(UCOr 

PSPaNNISE PRESSURE TaBLE SIZE EXCEEDED ' , FLOaT #NT A85C ) ) 
NTaBSC a NTaBSC - I 
00 379 I s J, NTABSC 
579 ThETsC(I) 8 0AHC0S(THET5C(U) 

GO TO 200 


NeTLOADS--«COMPUTES 4 ND PRINTS OVERALL RESULTS, CL,CBP 
(PITCH MOMENT), cMRfROLL MOMENT), LEFT aNO 
right ROOT RENqING moments, LEFT aN^ 

RIGHT LIFTS, LEFT aND RIGHT CENTERS OF PRESSURE, 
CDI, VORTEX ORaG FACTOR ( COI /CL**2/ ( P I * aR ) ) 

580 IF(CONV) write (UCO, 9012) 

REaO(UCI,6020J JJ,N0UT 
DO 399 L ■ 1 ,LL 
WRITe(U 6 f90l«) L 
C 

c the aspect Ratio in all the formulas worked out 
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c Has been the effective aspect Ratio, B**a/AREAf 
C but aR from GEOM file is «*BREF**2/AREA and BRATIO s 
2 * bREF/B, so the EFFErTIVE ASPEcT HaTIO is 
AREFF 8 AR/8RATIO**2 

also note that CBARBR 9 longitudinal reference length / 

BREF 50 CBARBa • CBARBR * BRATIO s 

longitudinal reference length / fB/2) 

STOHi « SRANlOaD C0EFFICIEnTS(BI,B 2.,,1 
CL*C/(2B) • 81 * SiNfTHETA) ♦ 8? • SIN(2*TH£TA) ♦ 

DO 3B1 K B I.KK 
STORl (K> a BNKI (li,K,L) 

581 IFfTESTNj) STOHlfK) b 8T0R1(K) ♦ ,5 * HNKir2,K,L) 

c lift coefficient 

CCL a CONI * STORl n) 

RR1TE{U6,90I6) CCL 
C 

«c induced ORaG 

c 

SUh a 0, 

K2 a KK A 1 
DO J82 K a 1,KK 
K2 a K2 • I 

582 SUM a SUM * FLOAT (K2 ) *STOR 1 (K2 ) **2 
COI 9 CONp A SUM 
HRITE(U6,9016) roi 

IF (STOsl f 1) ,LT, i ,F^2*5 UMi GO TO 365 
C 

c vortex orac Factor 

VORD a SUm/STORi f 1 }**2 
wRlTp(Ufe,9020) VDRD 
383 CONTINUE 
C 

c rolling moment 

CmR b C0N3*5T0RU2> / BRaTIO 
WRITEfU6,902l) CHR 
C 
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ROOT 8EN0ING MOMENTS 

CM8P B eCNOlNC MOMENT OF RIGHT HALF 

CMBM ■ sending MOMENT OF LEFT MaLF 

both OUaNUTIES are defined In the STRUCTURES SENSE, 

I,E,, rOTH will r€ positive FOR AN ELL IPT JC aLLY. 

loaded wiNr, with positive lift. 

SUM s 0, 

K B KK A 1 
DO J92 Kl » l,KK 
K s K * 1 

IF(MOD(R>2).EO,0) go to 592 

SACkwARDS summation USED FOR ACCURACY 
SUM s SUM + FL0AT({M0D(K*3»«)*1))/FL0AT(I<<‘*2-«) * 
ISTORI (K) 

392 continue 

CMPP * AREFF * (SUM+PIRASTORI (2) )/BRaTI 0 
CM6M « AREFF * (SUM^PlRASTORl (2 ) ) /SRaTTO 
WRlTE(U6,902a) CMRP 
WRITE(U6f9026) CHBM 


LIFT ON RIGHT 4NO LEFT SIDES 
SUM s 0, 

K ■ KK A i 
DO 393 Kl a 1,KK. 

K s K • I 

IF(M0D{K,2J ,NE,0) GO TO 39J 
SUM 9 SUM ♦ FL0AT(K*(M00rK+2, «),!))/ 
lFL0AT{K**2«n*ST0Hl (K) 

J9? CONTINUE 
C 

CLP a AREFFAfPI«*STORl{l)-SUH) 

ClH * AR£FFA(Pia*STORl f 1 )aSUh) 
WRITfc-<Ufr»9026) CLP 
WRITE (U6, 9030) CLH 



o ri o r% r> 


C SEMUSPaN centers op pressure 

TEST3«CLP.NE,0, 

1F(TESTJ)CPP« CM0P/CLP 
IP( TfSTS)WRlTE(U6r9032) c^P 
Test5«CLM,NE,0. 
IP(TEST3)CPMa-CM(i|«/CLM 
IP f TE5T3)«RlTEfU(!.,905«)CPM 


pitching moment 

pitching moment Can not be integrated analytically for 
A general planporm, subroutine intgwt will re used 
FOR the integration, 

OUM ■ AREPF/CBARB2 
DO J97 J ■ 8,JJMAX 

Sum « 0, 

SUMl B o', 

DO 396 K • l,KK 

inoek ■ J*P 

INDEX ■ MODf 1N0EX»1 f I8UB) ♦ 1 

ifmnoex,gt « jjmax]) go To 39a 

SKThET • SThETAHNoEx) 

CO to 395 

39« index ■ index , JJMAXI 
SKTHET b •STHET4(IN0EX) 

395 CONTINUE 

c 

C SKTHET B 51NfK*THIiiTA) WHERE 

C theta b ARC0S{ETA(JJ J 

C 

SUM a SUM • SKTHfT * STORifK) 

OUMi a BNKI n#K,U*5KTHET 

IP(TESTNn DUMl a DUMJ ♦ »NK 1 f 3, K, L > *SKTMET 
IF(T£STN2) DUMl 9 DUMl • ,S*BNK I ( 3 ^ K» L ) *SKThE t 
SUMl a SUMl 4 DUMl 

396 CONTINUE 
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SUM I •SUMMATION (K* 1 , KK ) S IN (K*THt'T A ^ * 

(BIK ♦ * b2K) s -CL*C/2 b 

SUM! s summation (K«t f KK)SIN(K*TH(;;TA)* 

fBlK ♦ B2K • ,5*81K) « -2,*CM(Lt) *C/8 
ST0R2(J) b OUM * tXSlLI»»(J)*SUH(.,25*CORDlP(J)*SUMn 
397 CONTINUE 

THE INTE0«ATI0n for PITCHING MOMENT WlLl BE DONE WITH 

increasing numbers of integration points, it will 
start with href points and increase up to JJMAX 
but not exceed jj. 

THE default for JJ IS JJMaX 

NOUT 18 the output LEVEL FOR SUBROUTINE 1NTg«T, 

NOUT B 0 IS USUAL 
WRITF(U6,9036) 

IFiJJ.EO.O) JJ B JJMAX 
J B MREF 
TeST3*CCL.Nt,0, 

DO 398 N » I, I ft 

call INTc,RT(J, JJMaX, JJMAX, eTA,STHETA# ST 0R2,N0 LNDT,eTaD, 
ISTORl ,U6, NOUT, NLEFT, CMP) 

IF{TEST3) XCPb-CMP/CCL 
iFfTFSTS) WRlTE(Ub,9038) J,CMP,XCP 
iFf ,N0T,TEST3)WRne(U6,9038) J,CMP 
J 8 J ♦ J 4 t 
IFIJ.GT, JJMAX) GO TO 599 
IFIJ.GT.JJ) GO TO 399 
39fl CONTINUE 
399 CONTINUE 
GO TO 200 


PRE8SURES---COMPUTE8 AND PRINTS ThE LIFTING PRESSURES At 
THE CHOROwISE LOCATIONS DEFINED 5V THE XCP 

command and the SPaNWISE locations defined 
by the ecp command, 

aoo continue 

DO «26 L«l,LL 
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w«ITF(U 6,90U) t 
00 NS r I ,NTABSC 

FTAS tt 0C0S(THtTSC (NS) ) 

IF( 5L0a 0<U , AND, ETAS, LT.-l ,£•<!) 00 TO «^5 

if(asloaO(U.ano,etas,lt, ue*«) go to uzs 

ETAS 


A LINEaK interpolation will rc used to obtain the 

CMORD/(B/g) /iT THE CURRENT ETaS STATION 

(ETa) is In oescenoino order 


IF(ETAS*LT,ETA(2)w1«D"(*) 00 To ttOl 
J] a ^ 

go to «o 5 
aol continue 


00 Uoi J a 3 ,JJHa;k 
J 1 * »I 

IF(ETaS,GT,ETA(J)) go to «03 
402 continue 
403 CONTINUE 


411 


412 


4,1 S 
‘ix'f 


C0R0B2 a CORDIP(Ji) ♦ (COROIP( Jj ) -CORD IP ( J !• i ) ) / 
I tET 4 (jn-ETA(Jl«'l))*(£TAS*ETA(Jl)) 

00 411 N • 1,NmAX 
STORj ( n) • 0 , 

00 412 K a 1 ,KK 


DTHETA • OSlN(OFLOAT(K)*THETSC<NS)) 

00 4i2 N a 1 ,NmAx 

STORKN) a STORlfN) ♦ 8 NkI(N,K,L) * OTHpTA 
00 424 NC a l^NTABC 
THETA K THETCW(NC) 

DELCP a STORl(l) /TaN( ,5*TheTa ) 

IF (NMax,LT,2) GO TO 417 
00 4!*> N a 2,NmAx 

DELCP X DfcLCP ♦ STORl(N) * OSIN(DpLOAT f N-l) *THfr TA ) 


DELCP « eOPi • 0ELCP/C0RD82 
X a •OcOS{THETA) 

CHI a ,5 * ( I , ♦ X) 
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WftITK (U6,90a«) K,CHl,OELCP 
aga CONTINUE. 

<l25 CONTINUE 
age CONTINUE 
r,0 TO 200 


FL4P 

aan continue 

WRiTf:(U6,9o«6) 
GO TO 200 


CONTINUE*— Causes program to continue execution in batch 

MODE EVEN IF aN invalid COMMaNO IS ENCOUNT£«ED, 
aeo CONTiN a .TRUE, 

GO TO 200 


STOP— halts execution 
aao write(U6,9050) 


STOP 


NEW. .starts a new case, REaOS N£W SOLUTION FILE aNO 
(MAYBE) A NEW GEOMETRY FRE 
485 RETURN 

KNEW--STARTS A NEW CASE, READS A NEW SOLUTION FRE 
and (MAYBE) A NEW geometry FRE, RETAINS CURRENT 

stations for pressure and spanloaos 
aoo return 
c 

500 continue 
C 

PWEIGHTS command PRINTS THE CURRENT WEIGHTS 

WRITE(UCO»8ott5) 

DO 505 Lai ,U 
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505 


w«lTC(UCO,^05n L# (wClGHT(I,L)f I«lfNC0lM) 

GO TO ?00 

— — * 

c 

520 continue 

I«.I«TS$ command allows user to give tss commands to operating system 

TSSCOM is an array In which to store the command and is not needed 

,,,,,0TMERWISE, 

read CUCI#R052) TSSCOM 
CALI OBEY(80f TSSCOM) 

GO TO ?00 


eOOO FORMAT!/, » MACH a •,F8,«/,* (B/2)/8hEF « », 

lFS,«/®» CBAR/BREF • t^PS a/, 

?» aspect ratio » «,F8,a/) 

8010 FORMAT!* * I ) 

8015 FORMAT!* ENTER COMBINATION CODE* ) 

8020 P0RMATM615) 

8025 FORMAT!* COMBINATION »,IJ/) 

6030 FORMAT!* ENTER WEIGHTS OF FIRST *,l3,l CASES* ) 

8035 FORMATfSplO.O) 

8000 FORMAT! I Inter weights of last *,I3,* CASES’ ) 

80R5 FORMATtlHl//* WEIGHTS'/) 

6060 F0RMAT!RA3) 

8065 FORMAT!* INVALID COMMAND *,«A3/) 

80T0 FORMAT!* ENTER NTYPE* ) 

80T5 FORMAT!* NTYP| INVALID*/) 

8060 FORMAT!* ENTER NO* ) 

6085 FORMAT!* ENTER NSTa* ) 

8087 FORMAT!* ENTER ETMIN, |TMaX , OET A * ) 

60RO FORMAT!* ENTER TaBlE OF ETaS— I PER LlNf *, 

1 ’ENDING WITH VAlUIE GREATER THAN i* ) 

8091 FORHAT!* enter TabLE OF X VALUES' ) 



6092 FORMtT(f ENTER TaBLE OF CHf VaLUES* ) 

8095 FOHMaT(« enter XMI N, XM a X, DX » ) 

9000 format t ETA <,7X1 ♦ CL*C/2B * , ‘CL,»#5X# 

I 'CL*C/CCe*CAVC #2X, »2Cm< 1/«)C/R| r2X, iC,P,(l/«) • ) 

9005 format ( iXf 4(F10,6r IX) >3Xr2fF10,6r IX) ) 

9012 FORMAT(t ENTER JJ iNO HOUT» ) 

90l« formats/, I combination ',13/ 

I I f) 

9016 F0RMAT(» Ci.‘i 17X#*« »#Fll.6) 

9018 FORMAT(f INDUCED DRAC » UFll.B) 

9020 FORMATtt VORTEX DRAG FACTOR a i,Fll,h) 

9021 FORMaT(» HOLUNG MOMfNT H/ (Q*8*2*BR£F ) • i 1 3X, IHa , 

IF12.6) 

9024 F0RMaT(i right ROOT rENOINC MOMENT HrR/ (Q*S*BMEF ) a «, 

1F11.6) 

9026 FORMAT(t LEFT ROOT rENQING MQMeNT MbU/ rO*S*»REF ) a •, 

IFl 1 ,6) 

9028 FOHMAT(i lift ON ETA.GT.O a »,Fll,6) 

00 R030 FOHMaT(« lift on ETA.LTjO » »,F11,6) 

« 9032 fohmaT(i Center of pressure of right half# v/bref a •> 

IF11,6) 

903« FOHMaTT* center of pressure of left half, Y/BREF b «, 

lFll.6) 

9036 F0HMAT(> PiTCHiNp MOMENT ABOUT X a 0'« 

i» and center of pressure »/ 

1» Ji, 7X, *CM/(0*S*C8aR) » ,7X, 'XCC.Pj/CeARf) 

9038 fORHAT( 3X, I3,FlB.9iFl9,6) 

9042 FORmaTC/* eta a * , F I I , 6 , / , 6X , l X ' , 1 2X , • CHI i , 1 0 X , 

1 'OELTA-CP' ) 

9044 format ( lX,plO,6,Fl 3,6,F 18.6) 

9046 FOR.MaT(/,» this FEaTU«E HaS NOT BEEN PROGS AMMEO </ ) 

9050 FORMaT(/, ‘PROGRAM HALTED DUE TO STOP COMMAND'/) 

9051 eORMaK ■ OCOMBINATION f , I3/(6 f10,6) ) 

9052 FORmAT(20AU) 

end 

subroutine INTGRT (JJ, JJH, JJMaX, eta, STHETa, F, no, X^ Y,4, 
lOUTR, NLEFT, VaLU) 

cccecccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c c 



orioo ririonoririnnorionrioririoonn 


C program To iNTtGRAtE PcX)*8tN()<) PROm 0 TO PI dV A TRAPEZOIDAL RUiE 
C M00IPI£D TO allow for SPECIFIED DISCONTINUITIES, 

c the discontinuities may be discontinuities in slope or Value 
c IF any oiscontinuitv is a discontinuity in Value and the location 

c OF THE DIScONlUITY COINCIDES WITH AN INTEGRATION POINT, THE PROGRAM 

C assumes THaT The corresponding element of f is filled with the limit 
c average value, 
c 

THE program is written TO BE MaCHInE INDEPENDENT 


JJ ■ number of stations To use in the integration 

jjM • number of Values of p 

JJMAx • number of Values of eta and stheta 

• COS(THETA) at jjmax stations 

SThETA ■ SiNftHlETA) AT JJMAX STATIONS 

F t THE FUNCTION VALUES AT JJR STaTIONS 

NO B THE NUMBER OF DISCONTINUITIES 

X B A vector CONTAINING THE DISCONTINUITIES 

Y B THE X VECTOR AFTER REORDERING FROM HIGH TO LOW aNO 

eliminating Values out of range, duplicated values, 
AND Values within eps of any integration station 
OUTR B AN integer governing TH£ OUTPUT OF THIS PROGRAM’ * 

nleft b the number of y. Values, 

VAlU b the output value of the integral 


ALL the input variables aRE RETURNED 


THE X. Values are completely unrestricted, hut the corrections will 
NOT necessarily be done correctly in all cases’. 


cccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c..,, .specification statements 

double precision ETa, sum, pi, to , SThETa 
integer OUTR, out, W 
logical CONT, COnTI 


cccccccccccccccccccccccccccc 


rt.r»dnr>nnr»r*r»oririnnnnrf<^r»ri<nr%r»-r»riori 



OIhENSION ETaC'I'JhAX), STHCTacJJmAX), F(JJR), X(nO), Y(NDJ 
nIMENSlON ETf«)» FUN(a) 

EOUIVaUENCE fYTEMP,VOUM) 

c 

DATA PI /S'.iai59265J58«»793 DO/ 

data EPS /t»E«5/ 

c 

JJ«, and jjmax must line up propehly, ip They dO not the 
c t » f t »**’^OCRaM KilUL write a MESSAtiE TO THaT EFFECT A*YD HaI,T EXECUTION 
c 

IF (MOO(JJMaX+I, JJR+J) ,N£. 0) GO TO 910 
IF(M0D(JJR ♦!» JJ fi) ,N|, 0) GO TO 9jO 

JRl «(JJRtn/JJl 
J«2 «(JJMAX*1)/JJ| 
delta 8 Pl/nFLOAT(JJl) 

SUMrO 
OUTb outr 
00 C 

'"C.,.,,THE 50 LOOP does the UNCORRECTEO integration 

c 

00 so 

Jt>J*JRl 

J2aJ*JR? 

SUm»SUM*SThfTA(J 2) *F f Jl) 

50 CONTINUE 

VALU a OElTA*SUM 
C 

there discontinuities, transfer to 200 

c 

ifcnd.gt.O) go To 200 
60 NlEFT«0 

IF( 0 UT,LT, 1 ) return 
wRITf (W^nVALU » JJ 

IF(0UT ,lt, 2) Return 

IFfND ,EQ, OT CO To 7o 

write (W»2) 

wRITE(W,3) X 

70 IF(0UT.LT,3) return 



WRITF (W,«) 

D075J®tfJJ 

75 WRITE (W,5) FTA(J?>, F(Jl) 

RETURN 

200 NLEeT a ND 
C 

C,,,,, PRIOR TO correcting FOR DISCONTINUITIES IN SLOPE aNO VaLUE THE 

EXTR4NEOU8 POINTS 4 RE eliminated from consideration, a transfer 

c back to so wji,l occur if There aRE no real Discontinuities. 

c 

c 

C,,,,,ARRaNGIN 6 K in y in descending order 

c 

do 210 Nal.MO 
210 V(N)a X(N) 

DO 220 Nal.NO 
DO 220 jaN.NO 
IF{Y(N),GE.Y(J))60 to 220 
VTEMPaY(J) 

V(J)aY(N) 

VfN)*YTEMP 
220 CONTINUE 

IF fNLEFT.EO.lJGO Tn JOO 
LMIN*1 
LMAxbNO.I 
C 

c repeated. Values are eliminated next 

c 

230 DO 200 LbLMIN^LHAX 

IF ( aBS< Y(L)»va*i ) ) .LT.EPSjGO to 25o 
200 continue 

GO TO 2feO 

250 Call crunch fY,L,.LMAXj 
IF(L.GT,LMaX) gO TO 260 
LMIN.L 
GO TO 23o 



260 NlCFT a U^AX+l 
300 CONTINUE 
C 

TOO shaU V. values A»e eliminateo next 

c 

YDUNa EPS*|t 
DO 3l0 Usl.NttFT 
IF( Y(L) .LT.VnUM) GO TO 320 
310 continue 
GO TO 330 
320 NLEFTat*! 

1F(NLEFT,E0,0) go to 60 
330 IF(Y(NLEFT) ,GT,»Y 0 UM) GO TO 60 

aoo IF(Nleft, eo. U Go to soo 
LMAXaNLt^T-l 
YbUMs.YDUM 
C 

C TOO LARGE Y-VALUeS ARE ELIMINATED NEXT 

« C 

'^iiio lF(Y(l),LE,YnUM) GO TO <U0 
CALL CRUNCH(Y, 1 ,LMaX) 

60 To «10 
NLEFT a LMAX + 1 
500 LMlN#i 
C 

C,, ,Y,VaLUES nearly coincident with integration stations are eliminated next 

c 

00 530 J«JR2,UJMAX, JR2 
ES a ETA(J) 

OO 510 LaLMIN, NLfFT 
IF(ABS(Y(l)-ES),Lt. EPS) GO TO 520 
510 CONTINUE 
GO TO 530 

520 CALL CRUNCH fYfLi NLEFT) 

IF (NLEFT ,eO*0)GO TO 60 
IFfL.GT,NLFFT )60 To 5«0 
LMJN a L 
530 CONTINUE 
C 



c at this point all the ext«aneou8 discontinuities Have been 

c.,»,,elihinated and There aRE still some remaining, the remaining dis- 

C..., .continuities a«E stored in (Y), 

c 

SRO KOUNT • 0 

CONT • , FALSE, 

VALUl»VALU 

IF(JJ,LT, «) GO TO 800 
C 

c..,..aT This POINT THE program hill Make the necessary corrections 

C it WILL NOT be able to DO THIS RIGHT UNDER CERTaIN CONOITIOnsI 

c...,,iE There are discontinuities in t^e first or last two intervals, 

c OR more than one discontinuity in an interval, oh discontinuities 

C,,.,,1N adjacent intervals, the corrections will not fit DONE RIGHT. 

Ef n)«ETA(JR2) 

ETt2)«ETA(2AJR2) 

ET(3)«ETA(J*JR2) 

CORRsO 

CONTI*, FALSE, 

JMAXbJJ.3 

LMIN.I 

DO TOO JaI,JMAX 
INOEX»f JY3)*JR? 

ET(«)«ETA(INnEX) 

DO 600 L» LMIN,NLEFT 
1F(Y(L) ,0T,FTf2:i) GO TO 600 
iFfY(L) ,CT,FTr3)) GO TO 6i0 
CoNTlw, FALSE, 

GO TO 690 
600 CONTINUE 
CO TO 79ft 
6t0 INOEx»ft 
JP3BJ+3 

00 6^0 KaJ,JP3 
K 1«K* JRI 
K?*K* JR? 

INoEXcINdEXaI 

620 FUNf INpEy)eF(Kn*5THETA(K2) 



TO c OARCOS(OBLE(YfU)) 

02» TO/OeLTA-OFLOAT(j4.1 } 

05 ■ 1, .02 

CORR B CORR ♦ D2**2*fFUN(2)-FUN( 1 > ) ♦ n3**2*(FUN(3)»FuN(«) ) 

1 ♦ (2. *02 • 1 ,)*rFUN(2)-FUN(5)) 

COnT • CONT^OR, CONTI 

COnTU .true, 

KOUNTsKOUNT+1 

lFtKaUNT,EO,NUFFT)GO TO 790 
690 continue 

ET(n»ETf2) 

ET(2)bET(3) 

ET( JJsET ftt) 

LMlNsU 

700 CONTINUE 

79o VAtUaVALU ♦CORRaOEL T A/2 , 

C... 

c 

C,,,,,TH1S Is The end of the computation, the remainder of the RROCRaM 

c,,,..18 OUTPUT, 

c 

800 IF (OUT ,LT, 0) RETURN 

IF (FOUNT, NE,NLEFT)G0 TO 0lo 
IF(cONT) GO TO 810 
lF(JJ,LT,a) GO TO elO 
IF (OUT.GT.O) go TO 820 
return 

810 OUT aMAXo(OUT»2) 

820 WRITf(W,63 VALUI.VaLu# JJ 
If (OU7.LT.2) return 
IF (FOUNT, fO.NUEFTj GO TO 830 
fount e NLfFT » FOUNT 
WRITp(wr7) FOUNT 
830 IF (.NOT, CONT) GO TO 8<i0 
WRITF(W»8) 

880 IF( JJ.GT,3)G0 TO 8S0 
RRITF(W,9) 

850 WRITE(W»10) 

WRITf(W,3) F 



WR1TE(w,U> 

W»ITP(W»5J (V(N),Nal ,N lEFT) 

RO TO 70 

910 WRITE fW, i2) JJ, JJR, JJMaX 

STOP 

1 format fgSHoVALUE OF TH£ INTEGRlt IS ,Fl5,7 iOX, 

i«0H THE number 0|P INTEGRATION POINTS USED a ,fS) 

? format JSOHONONE OF THE X^VaLUES W£»e CONSIDERED DISCONTINUITIES /. 

117H THE X VECTOR IS ) 

S format f 6F20,6 ) 

a format (loMo eta ^ lox, ihfj 

5 format ( F12.6, pn.fe) 

6 format t 

IttlHOTHE UNCORRECTEO VALUE Or THE INTEGRAL » I 

aaiM The corrected value Or the integral s .pis, 7/ 

SaiH the number Op integration points has a ,17 5 

7 format (2SM0-**WaRNING*»«THERE weRF , IS, 

152m discontinuities Which could not me corrected por, i’ 

y, 8 format ( eOHO«*«MARNlNGaw*THERE WERE DISCONTINUITIES IN CONTIGUOUS, 

® 1 INTEGRATION INTERVALS > 

9 format ( lOOHO— -W aRNING-**NO DISCONTINUITIES ARE CORRECTED FOR WHE. 

IN There are less Than a intec^ration points i 

10 format (26hoorici:nal discontinuities ) 

It format (25horeouceo discontinuities ) 

12 FORMAT fSSHoJj, JJR, a^O JJMaX aRE NOT PROPERLY REUaTEO /. 

19h0 JJ b , |7/ * 

19h JJR a , 17/ 

19h JJMav a , J7/ ■ 

I jaHoeXEcuTiON terminated in intgrt ) 

End 

subroutine crunch (Y, L» LMaX) 
dimension Yn) 

DO 10 k « Lf LMAj( 

10 Y(K) aYfK*l) 

LMAX • LMAx • I 

RFTURn 

End 

subroutine 8T0P2fN,ME8ACE» VAL) 
dimension MESAGEf2o) 






»RITC(N,i> «C8AGC,VAI 
STOP IS 

t FOPMAT(///i ***** ua.0*«i..» *****’/».. V*t.**fJPFi5. 71 . 



End 




